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ABSTRACT 

Context. Edge-on galaxies can offer important insights in galaxy evolution as they are the only systems where the dis¬ 
tribution of the different components can be studied both radially and vertically. The HEROES project was designed to 
investigate the interplay between the gas, dust, stars and dark matter (DM) in a sample of 7 massive edge-on spiral galaxies. 
Aims. In this second HEROES paper we present an analysis of the atomic gas content of 6 out of 7 galaxies in our sample. 
The remaining galaxy was recently analysed according to the same strategy. The primary aim of this work is to constrain 
the surface density distribution, the rotation curve and the geometry of the gas disks in a homogeneous way. In addition 
we identify peculiar features and signs of recent interactions. 

Methods. We construct detailed tilted-ring models of the atomic gas disks based on new GMRT 21-cm observations of 
NGC973 and UGC4277 and re-reduced archival Hi data of NGC5907, NGC5529, IC2531 and NGC4217. Potential 
degeneracies between different models are resolved by requiring a good agreement with the data in various representations 
of the data cubes. 

Results. From our modelling we find that all but one galaxy are warped along the major axis. In addition, we identify 
warps along the line of sight in three galaxies. A flaring gas layer is required to reproduce the data only for one galaxy, but 
(moderate) flares cannot be mled for the other galaxies either. A coplanar ring-like structure is detected outside the main 
disk of NGC4217, which we suggest could be the remnant of a recent minor merger event. We also find evidence for a 
radial inflow of 15 ± 5 km s“^ in the disk of NGC 5529, which might be related to the ongoing interaction with two nearby 
companions. For NGC 5907, the extended, asymmetric and strongly warped outer regions of the Hi disk also suggest a 
recent interaction. In contrast, the inner disks of these three galaxies (NGC 4217, NGC 5529 and NGC 5907) show regular 
behaviour and seem largely unaffected by the interactions. Our models further support earlier claims of prominent spiral 
arms in the disks of IC 2531 and NGC 5529. Finally, we detect a dwarf companion galaxy at a projected distance of 36 kpc 
from the center of NGC 973. 

Ksy words, galaxies: ISM - galaxies: structure - galaxies: kinematics and dynamics - galaxies: interactions - galaxies: 
individual 


1. Introduction 


Galaxies are dynamic, continually evolving systems in 
which the various components never cease to interact and 
influence each other’s evolution. Primordial atomic gas 
condenses into clouds of molecular gas, which then cool 
and collapse to form stars. At very low metallicities, stars 


can even form directly from the atomic gas (e.g. Krumholz 
2012[). Metals are produced by nucleosynthesis in stel 


lar cores and are injected into the ISM through stellar 


* Herschel is an ESA space observatory with science instru¬ 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 


winds and supernova explosions. These metals then either 
mix with the interstellar gas or condense to form small 


dust grains through various processes (see e.g. 

Dwek & 

Scaloj 1980 

[Kozasa et al. [19911 [Dwek||1998[ [Z 

lukovska 

& Henning 

2013 [Zhukovska 2014|l. On the other hand. 


dust grains are also continuously destroyed in the ISM, 
for example by supernova shock waves or processing by 
UV or X-ray photons, hot electrons or cosmic rays (e.g. 
JonesI 20041 [Micelotta et aLjpOIOl [Jones & Nuth|[^II[ 


Bocchio et al.| 20I2i. In their turn, dust grains stimulate 
star formation by catalysing the formation of H 2 molecules 
on their surface and shielding the molecular gas from the 
dissociating stellar UV radiation (Gould & Salpeter|1963 
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Hollenbach & Salpeter 1971[ Cazaux & Tielens 2002[ ). 
Dust grains are also efficient coolants of hot gas through 
inelastic collisions with thermal particles (e.g. Dalgarno & 
McCrayl 1972 Ostriker & Silk|T97^ . In addition to this 


‘closed-box’ evolution, galaxies also accrete external gas 
through merger events or gas infall from the intergalactic 
medium (e.g. ISancisi et al.|200^|Ocvirk et al.|2008j|Keres j 


et al.||2009b| Di Teodoro & Fraternali|2014| l, while super¬ 

novae and active galactic nuclei (AGN) can expel material 
from a galaxy (e.g. [Keres et al.|2009a[|Alatalo et al.|20k 


Lagos et al.||2013l [Sharma & Nath 20131 ). On top of this, 

according to the ACDM model, all these processes take 
place within a halo of dark matter (see e.g. Zwicky|1937 


|Trimble|19871 |Sahni|200^ , of which the nature and distri¬ 
bution, and the very existence, are still a matter of debate 


(e.g. |Milgrom|1983 de Blok|2010[ Feng|20 1 Oj |Famae}r& 


McGaugh|2012[ ). In order to uncover how galaxies form and 
evolve, it is therefore crucial to study all the components of 
a galaxy and the links between them, both locally and on 
global scales. 


Over the years many authors have investigated the con¬ 
nections between the different galaxy components. These 
studies, however, mostly targeted low to moderately in¬ 
clined galaxies, where the radial distribution of the com¬ 
ponents can be derived easily and where structural features 
such as spiral arms or bars are clearly visible. Edge-on 
galaxies can offer an important complementary contribu¬ 
tion to our understanding of galaxy evolution. Indeed, this 
is the only orientation where the vertical distribution of the 
different components can be studied. Moreover, due to line 
of sight (LOS) projection effects, the stars, gas and dust 
can be detected out to large galactocentric radii and low 
column densities. Another advantage of edge-on galaxies 
is that they allow the study of dust both in extinction and in 
emission. Indeed, in edge-on galaxies the dust often shows 
up as dark lanes in optical images. This opens up the possi¬ 
bility to also constrain the dust content by fitting detailed ra¬ 
diative transfer models to optical and NIR images (Xilouris 


etal.|1997 

19991 Bianchi|2007 [Baes et al.|20111|De Looze 

et al.j 2012 

De Geyter et al. 

2013| jSchechtman-Rook & 


Taking advantage of the unprecedented capabilities of 
the Herschel Space Observatory (|Pilbratt et al. 2010[), 


the HERschel Observations of Edge-on Spirals (HEROES) 
project was designed to map the distribution of the dust, 
gas, stars and dark matter in a sample of seven nearby, mas¬ 
sive edge-on spiral galaxies and study the links between 
these components. This will extend the recent study of the 
ISM of NGC 891 by [Hughes et al.| ( |2014| ) to a larger sam¬ 
ple of galaxies. New Herschel observations of the HEROES 
galaxies and a qualitative and quantitative analysis of them 
are presented in Paper I of the project (jVerstappen et al.j 
[25T31 ). 


In this second HEROES paper, we present a detailed 
analysis of the atomic gas content of the galaxies in our 
sample, based on a combination of new and archival in¬ 
terferometric Hi 21-cm observations. This analysis is cru¬ 
cial to the HEROES project for different reasons. Indeed, 
apart from the amount and distribution of the atomic gas 
in the disk of a galaxy, interferometric Hi observations 
yield a great deal of additional information. As the Hi gas 
is generally the most extended (visible) component, past 
or ongoing galaxy interactions are often most clearly or 
even solely visible in the atomic gas, for example through 


large-scale warps, a disrupted Hi morphology, tidal struc 
tures or the presence of companion galaxies in the same 
velocity range (e.g. Kamphuis & Brig^|1992| Yun et al. 


|1994t |Sancisi|ri999[ |Holwerda et al.||20l i| l. In the current 
work, NGC 5529 is a nice example of this. In addition to 
this morphological information, interferometric Hi obser¬ 
vations also allow to constrain the kinematics of a galaxy 
out to large radii, which in turn is necessary to derive the 
distribution of the dark matter. In order to unambiguously 
extract all the information that is available in the observa¬ 
tions, we perform a detailed tilted-ring modelling of the full 
data cube of each galaxy in our sample. The gas disk of 
each galaxy is represented as a set of concentric rings char¬ 
acterised by a number of geometric and kinematic param¬ 
eters, which are optimised to give the best possible agree¬ 
ment with the data cube. While time-consuming, this is the 
best strategy to resolve the degeneracies between different 
models of the atomic gas disk, which are especially present 
for an edge-on geometry. 

The Hi content in some of our target galaxies has al¬ 
ready been analysed in varying degrees of detail in the past 
(see Section |^. However, these studies often overlooked 
important information in the data cube or did not adapt their 
methods for an edge-on geometry. In order to obtain ho¬ 
mogeneous and accurate results for our entire sample, we 
re-analyse these galaxies along with the others. 

This work is the second step in our effort to obtain a 
detailed overview of the distribution and the properties of 
the different constituents of the galaxies in the HEROES 
sample. As a following step, we are gathering CO line ob¬ 
servations to also constrain the molecular gas component. 
Our team is currently also htting radiative transfer models 
to a set of optical and near-infrared (NIR) images of each 
galaxy to constrain the distribution of the dust and the stars 
(Mosenkov et al., in prep). Taking into account the distribu¬ 
tions of the gas and the stars, the rotation curves obtained 
in this work will then be modelled to derive the distribution 
of the dark matter (Allaert et al., in prep). Einally, bringing 
together our analyses of the various components, we will 
thoroughly investigate the connections between them, both 
radially and vertically. 

This paper is outlined as follows: in Section]^ we de¬ 
scribe the sample selection criteria, the observations and 
the data reduction. In Section[^we explain the general mod¬ 
elling strategy and our method to estimate the uncertainties. 
In Section]^ we present the results of our modelling and a 
short discussion of the main features for each galaxy. In 
Sectionj^a constant scale height is weighed against a flared 
gas layer in the light of theoretical expectations. The ef¬ 
fect on the gas disks of interactions with satellite galaxies is 
briefly discussed in Section]^ and, finally, our conclusions 
are presented in Section|7] 


2. Observations and data reduction 

2.1. The sample 

The HEROES project targets a sample of seven massive, 
edge-on spiral galaxies. A detailed description of the sam¬ 
ple selection is given in Paper I, so we restrict ourselves to 
a short summary here. One of the goals of the HEROES 
project is to map the dust emission in the EIR/submm with 
the Herschel telescope and compare this to the predic¬ 
tions from detailed radiative transfer models fitted to op- 
tical/NIR images. To ensure sufficient spatial resolution in 


























































































Table 1: Hi observations 


Galaxy 

Telescope 

Array 

Date 

Bandwidth 

(MHz) 

Number of 
channels 

Time on 
source (h) 

Reference 

NGC 973 

GMRT 


Aug 2011 

4.2 

256 

5.7 

1 

UGC 4277 

GMRT 


Jul 2012 

16.7 

512 

6.2 

1 

IC2531 

ATCA 

6.0A 

Apr 1996 

8.0 

512 

9.9 

2 

IC2531 

ATCA 

6.0B 

Sep 1996 

8.0 

512 

10.0 

3 

IC2531 

ATCA 

1.5D 

Mar 1997 

8.0 

512 

9.1 

3 

IC2531 

ATCA 

750A 

Jan 2002 

8.0 

512 

9.6 

4 

NGC 4217 

WSRT 

Not mentioned 

Aug 1990 

5.0 

63 

12.0 

5 

NGC 5529 

WSRT 

Maxi-Short 

May 2001 

10.0 

128 

11.9 

3 

NGC 5907 

VLA 

Modified C 

Aug 1997 

2 X 3.1 

2 X 32 

4.7 

6 


References. (1) This work; (2)|Bureau & Freeman||(|1997); (3)|Kregel et al.|(|2004J; (4)|0’Brien et al.|(|2010d); (5)|Verheijen & Sancisi| 
pOOT) ; (6) |Shang et aLjp^ . 


all the Herschel bands, the galaxies were selected to have 
an optical diameter of at least 4'. They also need to have a 
well-defined and regular dust-lane to allow radiative trans¬ 
fer modelling. The latter criterion automatically rules out 
all galaxies with an inclination that deviates more than a 
few degrees from edge-on. Fitting radiative transfer models 
to optical/NIR images of galaxies was already done in the 
past by Xilouris et al. ( |1997[ 1999[ ) and Bianchi ( |2007| l with 
their respective codes. Their combined samples therefore 
formed the starting point for the HEROES sample. After 
applying our selection criteria, eight galaxies remained, of 
which NGC891 was omitted due to overlap with another 
Herschel program. An overview of the general properties 
of the HEROES galaxies can be found in Table 1 of Paper 
I. 


Although part of the HEROES sample, NGC4013 will 
not be treated in this paper. The Hi content of this galaxy 
was recently extensively modelled by |Zschaechner & Rand| 
( |2015[ ), according to the same strategy as followed in this 
work. In summary, they find that the gas disk of NGC 4013 
has both a substantial edge-on warp and a LOS warp. In 
addition, the gas disk is flared, with the scale height in¬ 
creasing from 280 pc in the center to about 1 kpc at radii 
larger than 7 kpc. Einally they report a lag (decrease of ro¬ 
tational velocity with increasing scale height) that becomes 
more shallow with radius. No extraplanar Hi gas is found, 
despite the presence of extraplanar dust and diffuse ionised 
gas. 


2.2. Observations 

NGC 973 and UGC4277 were observed with the Giant 
Metrewave Radio Telescope (GMRT) in August 2011 and 
July 2012, respectively. The GMRT consists of 30 45-m an¬ 
tennas located at a fixed position, yielding only one possi¬ 
ble antenna configuration. A central area of approximately 
one square kilometer contains 14 antennas, with the short¬ 
est baselines measuring about 100 m. The other antennas 
are placed in a roughly Y-shaped configuration with three 
arms of about 14 km long, giving a maximum baseline of 
25 km. In L-band (1000-1500 MHz), each GMRT antenna 
provides two orthogonal, linearly polarised signals. Both 
galaxies were observed for about 8 hours, including the 
observation of standard flux and phase calibrator sources. 
NGC 973 was observed with a bandwidth of 4.2 MHz, di¬ 
vided into 256 channels. Eor UGC4277, a bandwidth of 
16.7 MHz divided into 512 channels was used. 


Eor the remaining galaxies discussed in this pa¬ 
per, we used literature data from the Very Large Array 
(VLA), Westerbork Synthesis Radio Telescope (WSRT) 
and Australia Telescope Compact Array (ATCA). A brief 
overview of these observations is given in Table[2 We refer 
the reader to the original papers for further details of the 
data. 


The data for NGC 5907 and NGC 5529 were obtained 
from the VLA and WSRT archive respectively. Eor IC 2531, 
the flagged and calibrated Hi data used in the series of pa¬ 
pers by O’Brien et al.| ( 2010d|b|a|cl l and in his own work 
( [Peters et al.|2013[ ) were kindly made available by S. Peters. 
These consist of a combination of data sets obtained be¬ 
tween 1992 and 2002 with different configurations of the 
ATCA. Upon inspection of the individual data sets, how¬ 
ever, it became evident that two of them were severely cor¬ 
rupted, so we did not include them in our final data cube. 
Part of these data were also used by jKregel et al.| ( |2004| l; 
Kregel & van der Kruit|((2004|l. Einally, M. Verheijen was 


kind to provide a fully reduced and deconvolved data cube 
of NGC 4217 from his WSRT observations ( [Verheijen &| 
[Sancisi 1200T] |. 


2.3. Data reduction 

The reduction of the GMRT and VLA data (NGC 973, 
UGC Mil and NGC 5907) wa s performed in the AIPS soft¬ 
ware package ( [Greisen[|2003[ l following the standard pro¬ 
cedure. To save time in handling the large data volumes 
at hand, the initial part of the reduction was performed on 
the ‘channel-0’ data set (i.e. the average of the central 75% 
of the frequency channels). This data set was visually in¬ 
spected and data corrupted by radio frequency interference, 
malfunctioning hardware, etc. were flagged. Next, the flux 
and phase calibration was performed for the channel-0 data. 
Calibration of VLA or GMRT data in AIPS is based on the 
observation of a flux calibrator at the start and the end of the 
observing run and periodic observations of a phase calibra¬ 
tor during the run. Eirst, the observed flux calibrator was 
used to set the absolute flux scale. Next, the antenna gain 
amplitude and phase solutions were determined for both 
calibrator sources and the absolute flux information from 
the primary calibrator was transferred to the secondary cal¬ 
ibrator. Einally, the antenna gain solutions from the calibra¬ 
tors were interpolated to obtain the antenna gains for the 
target source. After the initial reduction of the channel-0 
data, the antenna gains and bad data flags were copied to the 





















































Table 2: Parameters of the final data cubes used in the modelling 


Galaxy 

Synthesised beam FWHM 

Velocity resolution 

RMS noise 

mjy/beam to Kelvin 


(arcsec) 

(kpc) 

(kms"^) 

(mjy beam~^) 

(K) 

NGC 973 

24.5 X 15.3 

7.5 X 4.7 

6.9 

1.6 

1.62 

UGC 4277 

15.6 X 9.8 

5.8 X 3.6 

6.9 

1.0 

3.96 

IC2531 

16.0 X 9.2 

2.8 X 1.6 

6.6 

1.1 

4.12 

NGC 4217 

18.0 X 13.0 

1.7 X 1.2 

33.0 

1.0 

2.59 

NGC 5529 

23.6 X 14.0 

5.7 X 3.4 

16.5 

0.6 

1.83 

NGC 5907 

14.8 X 13.3 

1.2 X 1.1 

20.6 

0.3 

3.08 


spectral line data, which were finally bandpass calibrated. 
Typically several iterations of flagging and calibrating were 
needed to achieve a satisfactory result. 

The WSRT data of NGC5529 were reduced with the 
Miriad software package ( |Sault et al.|1995 i. Although the 
global strategy of the data reduction is the same as in AIPS 
(flagging, calibrating and imaging the data), the exact cali¬ 
bration procedure is somewhat different. The observations 
of NGC 5529 did not include any periodic observations of a 
phase calibrator. Only a flux calibrator source was observed 
at the start and the end of the observing run. This was used 
to set the absolute flux scale and obtain initial antenna gain 
solutions, as well as to determine the bandpass corrections 
to be applied. To take into account the strong variations of 
the antenna gain phases during the observing run, the qual¬ 
ity of the initial calibration was improved through several it¬ 
erations of self-calibration. This was performed on the con¬ 
tinuum emission to increase the signal-to-noise ratio and 
the derived antenna gains were subsequently copied to the 
spectral line data. 

In order to suppress potential Gibbs ringing in the spec¬ 
tra, a data set is usually Hanning smoothed. For NGC 4217 
this step was already performed online, i.e. during the ob¬ 
servations. For the remaining galaxies this was not the case 
and the data had to be Hanning smoothed manually during 
the data reduction. The inevitable downside of this smooth¬ 
ing, however, is that adjacent channels are no longer in¬ 
dependent. One usually accounts for this by dropping ev¬ 
ery second channel, thereby reducing the spectral resolu¬ 
tion of the final data cube by a factor of 2. The purpose 
of our analysis is to fit 3D models to the Hi data cubes of 
the galaxies in our sample. This requires a good constraint 
on their kinematics and hence a sufficiently high veloc¬ 
ity resolution, preferably below 20 km s“^. For NGC 5907 
and NGC 5529, Hanning smoothing would lead to a res¬ 
olution of 41 km s“^ and 33 km s“^ respectively. We 
therefore decided to perform Hanning smoothing only for 
NGC 973, UGC4277 and IC2531 and not for NGC 5529 


and NGC5907. For NGC4217 the resolution is also 33 km 
s“^, but unfortunately nothing could be done about this, 
since Hanning smoothing was already applied during the 
observations. 

After subtracting the continuum in the uv plane, the data 
were imaged using a robust parameter of 0. This provides 
a good compromise between a high resolution and a high 
signal-to-noise ratio. The maximum baseline of the GMRT 
is approximately 25 km, which corresponds to a maximum 
uv distance of 119 kA at a wavelength of 21 cm. Using 
this entire uv range to image the data results in a very nar¬ 
row synthesised beam and produces a data cube where the 
emission is effectively resolved out. To image the data of 
NGC 973, the uv range was therefore limited to 40 kA and 
a Gaussian taper with a width of 10 kA at the 30 % level 


was applied. For UGC 4277, a better result was obtained by 
limiting the uv range to 15 kA and not applying any taper. 

The deconvolution of the dirty beam was generally done 
using the CLEAN task in Miriad. This task offers a choice 
between the Hogbom ( Hogbom|1974[ ), Clark ( |Clark|1980 1 


and Steer (|Steer et ak 19841 Clean algorithms and by de¬ 


fault determines by itself which is the most appropriate al¬ 
gorithm to use. After an initial Clean of the data cube down 
to 3 times the rms noise (a), the data cube was smoothed to 
a beam size of 60" and 2a masks were defined in the result¬ 
ing cube. The original data cube was subsequently Cleaned 
down to Icr in the masked regions and the resulting CLEAN 
model was then used as input model for a final Cleaning of 
the original cube down to ha. Einally, the Clean compo¬ 
nents were convolved with a Gaussian beam fitted to the 
central lobe of the synthesised beam and restored to the 
residual image. A different strategy was used for the de- 
convolution of the cube of NGC 5907, the closest and by 
far the largest (in terms of angular size) galaxy in our sam¬ 
ple. Ordinary Clean algorithms typically represent a source 
by a collection of point sources, which is a poor approxi¬ 
mation for extended sources like NGC 5907. To overcome 
this problem we used the multi-scale Cleaning option in 
AIPS ( |Greisen et al.||2009 i, with an enhanced version by 
Bill Cotton and Fred Schwab ( |Schwab|1984| l of the Clark 
Clean algorithm. Multi-scale Cleaning essentially uses a set 
of circular Gaussians as model components instead of point 
sources. In the current case three of these Gaussian com¬ 
ponents, with widths of 40", 120" and 360", were used in 
combination with a point source component and each chan¬ 
nel map was Cleaned down to 2a. 

Finally, each data cube was corrected for the primary 
beam and transformed to barycentric radio velocities. An 
overview of the properties of the final data cubes is pre¬ 
sented in Table |2] 

To produce total Hi maps of our galaxies, we first 
smoothed each cube to a resolution of twice the original 
(Clean) beam size. Masks were then generally defined in 
the smoothed cubes by excluding all emission that does not 
exceed 2a in three or more consecutive channels. If the ve¬ 
locity resolution of a data cube is low, however, three con¬ 
secutive channels can already span a large velocity range. 
For NGC 5529, NGC 5907 and especially NGC 4217, the 
latter criterion was therefore somewhat relaxed. After ap¬ 
plying the masks to the original data cubes, their zeroth 
moment was calculated. 

Following standard practice, we assume that the Hi 21 
cm line is completely optically thin. This means that the 
observed 21 cm surface brightnesses can be translated di¬ 
rectly into atomic gas column densities. In fact, it has been 
known for a fairly long time that this assumption does not 
always hold in the ISM of a galaxy. Observations of the 
absorption of continuum background emission by clouds 























of atomic hydrogen in the Milky Way have shown that Hi 
gas can indeed become self-opaque at 21 cm (e.g. [Hagen &| 
|McClain||1954[ |Radhakrishnan||l9Mt [Gibson et al.||2005f l^ 

The most detailed study of Hi self-opacity to date was car¬ 
ried out by [Braun et al.j ( [2009[ ) and |Braun| ( [2012[ l. They 
analysed high-resolution Hi data cubes of M31, M33 and 
the LMC and suggested that these galaxies contain a popu¬ 
lation of dense atomic hydrogen clouds that are self-opaque 
in the 21 cm line. These structures are compact, with typi¬ 
cal sizes of the order of 100 pc. Correction factors for the 
total atomic gas mass of 1.30, 1.36 and 1.33 were found for 
M31, M33 and the LMC respectively. However, since self¬ 
opacity occurs on such small linear scales, it is very difficult 
to correct for this effect in a meaningful way in more distant 
galaxies (observed with lower physical resolution). Studies 
of the atomic gas content of galaxies, both for moderately 
inclined systems and edge-ons, therefore generally assume 
that the 21 cm line is optically thin. 

An important question to ask now is what the conse¬ 
quences of this assumption are in an edge-on galaxy, where 
each line of sight cuts through the entire plane of the disk. 
Fortunately, for self-absorption to be effective. Hi clouds 
must not only overlap in projected position, but also in ra¬ 
dial velocity. Due to the rotation of a galaxy, the latter is 
generally not the case along a line of sight and the prob¬ 
lem is in fact reduced to the level of the individual channel 
maps. Hence the effects of assuming an optically thin 21 
cm line for edge-on galaxies will not be much larger than 
for studies of low or moderately inclined galaxies. 


3. Modelling strategy 

3.1. General strategy 


To model our Hi data cubes, we made use of the tilted- 
ring fitting software TiRiFiC pdzsa et'ar]|2007| l. For each 
galaxy, the width of the rings making up the models was 
taken to be half the size of the minor beam axis, as a com¬ 
promise between angular resolution and efficiency. Initial 
guesses for the relevant parameters were derived as fol¬ 
lows. The inclination was taken over from the literature. 
The (central) position angle and systemic velocity were es¬ 
timated by eye from the observed data cube. With these val¬ 
ues as inputs, we then used the task rotcur from the GIPSY 
software package (van der Hulst et al. [I992| l to perform 
a tilted-ring modelling of the observed velocity field. This 
technique is actually not well-suited for an edge-on galaxy 
and yielded rough initial estimates of the central postion 
and rotation curve. The radial surface density profile was 
derived with the GIPSY task radial, which essentially ap¬ 
plies the method described by [Warmels ( [1988| ) and [Lucy 


( [1974[ ). An initial value of 0.3 kpc was used for the verti¬ 
cal scale height and finally the velocity dispersion was es¬ 
timated as the squared sum of the velocity resolution of the 
obse rvations and an in trinsic velocity dispersion of 10 km 
s“^ (Leroy et al. 2008 1 . For NGC 4217, the results from the 
analysis by Verheijen & Sancisi|p001[ l were used as initial 
estimates. All models use a sech^ distribution for the verti¬ 
cal gas distribution. This profile was originally devised for 
stellar disks ( [Spitzer|1942[l, but is als o often used to model 
gas disks (see e.g. Moster et al.|2010[ Gentile et al.|2013 i. 

These initial estimates were used as input to perform 
a first TiRiFiC fit of the data cube. A very basic geome¬ 
try was imposed here, consisting of a one-component disk 
with a constant scale height and inclination across the entire 


disk. The first fit was performed with the central position, 
the systemic velocity, the velocity dispersion and the posi¬ 
tion angle as free global parameters, and the rotational ve¬ 
locity free for each ring. For warped galaxies, the position 
angles were then adjusted to match the warp. Subsequently, 
the results were used as input in a second fit, this time leav¬ 
ing the inclination and vertical scale height free as global 
parameters and fitting the surface brightness for each ring. 
The resulting basic model was then further refined by in¬ 
vestigating different geometries and gradually adding extra 
features, each time optimising the column densities and ro¬ 
tational velocities (if necessary) to obtain the best possible 
fit to the data. The approaching and receding halves were 
modelled separately, with the (obvious) constraint that the 
central position, the systemic velocity and the velocity dis¬ 
persion should be the same for both sides. In the basic ge¬ 
ometry the inclination and scale height were also required 
to be the same on both sides. The detailed strategy followed 
for each galaxy depends on the specific features encoun¬ 
tered and is described in the following sections. 

An automated minimisation is included in the 
TiRiFiC code, which is performed after convolving the 
model cube with the observational beam. Although this 
greatly expedites the optimisation of the parameters, its re¬ 
sults should not be trusted blindly. With a constant noise 
across the maps, by construction this minimisation scheme 
will assign more weight to regions of higher surface bright¬ 
ness. However, low surface brightness features often im¬ 
pose the strongest constraints and are most important to 
discriminate between models. Also, TiRiFiC does not take 
into account whether a model is realistic or not. The the¬ 
oretically best fitting model will often contain unphysical 
values and unrealistically strong radial variations of the pa¬ 
rameters. Therefore, the geometry of the models was gen¬ 
erally imposed by hand, and the parameter values obtained 
from a TiRiFiC fit were further refined by manual adjust¬ 
ment. 

The edge-on orientation of the galaxies in this work 
brings along a number of potential degeneracies between 
different models. A prime example of this is the vertical 
thickness of the gas disks, as seen in projection on the plane 
of the sky. This can be modelled by adjusting either the ver¬ 
tical scale height or the inclination of the disk, or a combi¬ 
nation of both. A second important example is the inter¬ 
pretation of the column density peaks observed in a zeroth 
moment map, where each line of sight contains the sum 
of multiple rings at different galactocentric radii. In order 
to resolve these degeneracies, it is important to use all the 
information provided by the observations and to compare 
the models to the data in different projections of the data 
cube. Scale height and inclination effects, for example, can 
be disentangled by studying the shape of the emission con¬ 
tours in the different channel maps. Increasing the scale 
height simply makes the emission thicker in all the channel 
maps. On the other hand, decreasing the inclination moves 
the emission away from the major axis and gives it a V- 
like shape in the channel maps (see e.g. Fig. [plater in this 
work for an example of this). The strength of this effect also 
strongly varies within the data cube, going from hardly vis¬ 
ible at high line of sight velocities to clearly visible close 
to the systemic velocity. Density peaks in a total Hi map 
can be understood by inspecting the associated major axis 
position-velocity diagram. The addition of the kinematical 
information allows to assign the correct surface density to 










































NGC973 


each ring and potentially identify additional features that do 
not belong to the axisymmetric disk (see e.g. Figs. and 
[plater in this work). 


3.2. Uncertainties 


Following standard practice, the modelling strategy that we 
use in this work is a combination of minimisation and 


the personal judgement of the modeller (e.g. Zschaechner 


et al.|2012| Gentile et al.|2013j|Kamphuis et al.|2013| l. As 

a consequence it is not possible to determine the errors on 
the parameters of the final models in a statistical way. We 
therefore estimated the uncertainty on each parameter by 
keeping all the other parameters fixed to their best fit val¬ 
ues and varying the parameter in question until the model 
clearly deviated from the observed cube. For each param¬ 
eter, there are usually certain regions of the disk and/or 
certain projections of the data that are especially sensitive 
to variations of that parameter. We focussed on these re¬ 
gions and projections to assess the deviation of the model 
from the observed cube. Errors on the central position and 
the position angle were determined based on the moment-0 
maps. For the inclination and scale height we focussed both 
on the moment-0 maps and on the inner and outer channel 
maps, respectively. For the systemic velocity we focussed 
on the central part of the major axis position-velocity (XV) 
diagrams and the uncertainties on the velocity dispersion 
were determined from the spacing between the contours on 
the terminal edges of the XV-diagrams. Finally, errors on 
the rotation curve were determined as the maximum of two 
values: half of the difference between the approaching and 
the receding side, and the channel width of the data cube 
divided by the sine of the inclination. 



RA (J2000) 


4. Modelling results 

4.1. NGC973 

NGC973 (UGC2048) is located at a distance of 63.5 
Mpcj^ It is classified as an Sb or She galaxy and has a 
boxy/peanut-shaped bulge ( |Lutticke et al.|200^ classify it as 
a ‘bulge which is on one side boxy and on the other peanut¬ 
shaped’). After subtracting an axisymmetric model from 
their I-band image of NGC 973, [Patsis & Xilouris| ( |2006| l 
find an X-shaped pattern in the central regions which they 
interpret as evidence for the presence of a bar. So far, the 
structure and kinematics of the atomic gas disk of NGC 973 
have never been studied. 


Fig. 1: Top: False colour V-band image of NGC 973 based 
on observations made with the Italian Telescopio Nazionale 
Galileo (TNG) operated on the island of La Palma by the 
Fundacin Galileo Galilei of the INAF (Istituto Nazionale 
di Astrofisica) at the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astrofisica de Canarias. 
The data were cleaned and combined using standard data 
reduction techniques. Bottom: Hi contours overlaid on the 
same V-band image. The contours start at 1.76 x 10^° atoms 
cm“^ and increase as 3, 6, 9, 12 and 15 times this value. 
The Hi beam is shown in the bottom-left comer. North is 
up. East is to the left. 


4.1.1. The data 

Our GMRT data of NGC 973 have a total integration time of 
5.7 hours, which is rather limited for a galaxy at a distance 
of 63.5 Mpc. Consequently the signal-to-noise ratio in the 
channel maps of our data cube is generally on the low side. 
As a result, the shape of the observed emission is rather ir¬ 
regular and sometimes varies significantly from one chan¬ 
nel map to the other. Unfortunately not much can be done 
to improve this. We therefore do not expect our models to 
nicely reproduce the shape of the emission in each channel 
map. 


' The distances we use here are taken from NED, as explained 
in section 2.1 of Paper I. 


The zeroth moment map of the Hi disk of NGC 973 is 
compared to a V-band image in Eig. The NE side of the 
disk is the approaching half, the SW side is the receding 
half. Both the stellar and the Hi disk are remarkably flat, 
even in the outer regions, but some asymmetry is visible at 
high contour levels in the total Hi map. This is also vis¬ 
ible in Eig. 1^ Indeed, on the approaching side only one 
density peak is observed (at a major axis offset of +\' in 
Eig.|g. This peak is mirrored on the receding side, at an 
offset of -I', but additionally a second, more central peak is 
also present on this side. The maximum (integrated) surface 
density is also significantly higher on the receding side. 

We further report the detection of a small satellite 
galaxy slightly above the NE edge of the disk, at a dis- 




















4.1.2. Models 
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Fig. 2: Total Hi maps of the various models discussed here 
as compared to the observed total Hi map of NGC 973. The 
B model is our best model. Contour levels are the same is 
in Fig. 


tance of 36 kpc from the center of NGC 973 (as determined 
from our models). We estimate its central position as RA 
02^34™23.9®, Dec 32°32'5.0" and find no optical counter¬ 
part. The satellite appears in 9 consecutive (independent) 
channels in the data cube, between 4517 and 4572 km s“^, 
although often only at the 3(T level, and has a total Hi flux of 
0.36 Jy km s“^. This translates to a Hi mass of 3.39 x 10® 
Mq. 


For the disk of NGC 973 we find an integrated Hi flux 
of 9.32 Jy km s“^, correspo nding to a total atomic hydro- 


gen mass of 8.87 x 10® M 0 . [Huchtmeier & RichtCT 
report a sing le dish flux of 8.^ Jy km s~^ while|Springob 


et al. (2005 1 give a correctecQ single dish flux of 13.9 Jy 


km s“^. Although these flux measurements are of the same 
order, they show a relatively large spread, which is caused 
by a combination of different effects. The dominant factor, 
in general, is the uncertainty in the flux calibration of the 
different observations, which can cause differences of up 
to a few ten percent. Secondly, given the limited observ¬ 
ing time, a significant reservoir of low column density gas 
might still be hidden in the noise in our channel maps. In 
addition, single dish observations often have a half power 
beam width (HPBW) that is smaller than the HI disk of the 
galaxy, causing them to miss some of the emission from 
the gas. To correct for this effect, Springob et al. ( |2005| l 
assume a smooth Gaussian distribution for the HI, scaled 
according to the optical radius. Given the variety of mor¬ 
phologies shown by the atomic gas disks of galaxies, such 
a general correction can also be rather inaccurate for a spe¬ 
cific galaxy. 


For the modelling of the data cube we focus only on the Hi 
disk of NGC 973. The satellite galaxy is not included in the 
models. The best fitting model in the basic geometry (which 
we will refer to as the basic or B model) has a constant incli¬ 
nation of 88° and scale height of 0.2 kpc, although the latter 
is very ill constrained, with values up to 1.1 kpc giving only 
marginally different models. A single, constant position an¬ 
gle of 228° was found to be sufficient to reproduce the ori¬ 
entation of the disk on the plane of the sky. If we compare 
the zeroth moment map of this model to the observations 
(Fig. 13 we see that the agreement is generally very good, 
apart from the observed central density peak on the reced¬ 
ing side which is not present in the model. The main trend 
in the channel maps is also well reproduced by the model 
(Fig. |3 I, although (as expected) it often does not match the 
exact shape of the observed emission. Additionally, at inter¬ 
mediate rotational velocities on the receding side the model 
often underestimates the surface densities and hence lacks 
the high surface brightness contours. This effect can, for ex¬ 
ample, be seen in the channel maps at 4565,4703,4840 and 
4909 km s“^ and is also visible in the major axis position- 
velocity (XV) diagram (Fig. |3- The latter is well modelled 
at the terminal sides (i.e. at the maximal velocities) and in 
the very central region, but lacks high surface brightness 
emission at intermediate velocities, especially on the reced¬ 
ing side. In a regularly rotating, axisymmetric disk that is 
seen edge-on, the emission that is seen at these intermedi¬ 
ate velocities is simply the projection along the major axis 
of the emission at the terminal velocities. Hence, the only 
ways to significantly enhance this emission are either to in¬ 
crease the inclination of the disk or to increase the surface 
brightness in part or all of the rings making up the model. 

The first option was tested by constructing a second ba¬ 
sic model, the B2 model, with a constant inclination of 90°. 
The corresponding optimal scale height is 0.4 kpc. With the 
surface brightnesses optimised to match the terminal sides 
of the XV-diagram, this model offers little improvement on 
the receding side (Fig.|3- On the approaching side the prob¬ 
lem is almost solved, but at the cost of now overestimating 
the emission close to the systemic velocity. 

In a second attempt to improve our basic model we 
manually scaled up the surface brightness distribution on 
the receding half of the B2 model until we could repro¬ 
duce the high contours at intermediate velocities in Fig. |3 
With an increase of 35 percent {B3 model) this is almost 
achieved, but the XV-diagram is now clearly too wide on 
the receding side. This effect can also be seen in the chan¬ 
nel maps at 4909, 4977 and 5046 km s“^ in Fig. where 
the left and right edges of the model extend well beyond 
those of the observed emission. In addition, the emission 
at the terminal and systemic velocities in the XV-diagram 
is strongly overestimated (note that the B3 model includes 
a 6cr contour at the systemic velocity in Fig. |3 while the 
observed XV-diagram does not). In the total Hi map (Fig. 
13 the B3 model overestimates the peak density on the re¬ 
ceding side and shows one big peak instead of two separate 
peaks. 

From the B2 and B3 models it is clear that the Hi disk 
of NGC 973 cannot simply be modelled as an axisymmet¬ 
ric disk. Either we underestimate the amount of gas moving 


^ Corrected for pointing offsets and source extent, not Hi self¬ 
absorption. 





































Data 


B: Basic 


B2: IncI = 90° 


B3: Extra gas at rec. side 
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Fig. 3: Representative channel maps from the observed data cube of NGC 973 and the various models. Contour levels are 
-2.3, 2.3 (l.Str) and 6.2 mJy beam~^. The observations are shown as black contours, with the negative contours as dashed 
grey. The green contours show the final (B) model. Other models are shown as red contours. The line of sight velocity is 
indicated in km s“^ in the bottom left comer of the first panel of each row. The systemic velocity is 4770 ± 8 km s“^. 
Note that the horizontal axis is slightly shifted between the three blocks composing the figure. 


at intermediate velocities or we overestimate the amount of 
gas moving at the terminal velocities. The reason for this is 
most probably that the disk of NGC 973 is far from smooth 
and has a clumpy structure with high-density regions al¬ 
ternated by regions with a lower density. Alternatively, a 
projected spiral arm might also cause this behaviour. In any 
case, although TiRiFiC allows to model different areas of a 
gas disk individually and to include spiral arms, the prob¬ 
lem at hand, in an edge-on disk, is much too degenerate 
to attempt such a thing. In addition, apart from the arbi¬ 
trary jumps of the observed emission in the channel maps. 


we do not see a clear trend indicating that the shape of the 
model emission is systematically wrong. Hence, we stick to 
the basic geometry and use the B model as our final model 
for NGC 973. We stress that this uncertainty concerning the 
exact morphology of the disk does not strongly affect our 
results and is primarily restricted to the geometric parame¬ 
ters, i.e. the inclination and the scale height. Slight adjust¬ 
ments of the surface density profile can be necessary when 
adding a flare or a line of sight warp to a model, but these 
are typically only of the order of 10% or less. In principle 
a LOS warp also affects the rotational velocities, but for 



































B: Basic B2: IncI = 90° B3: Extra gas at rec. side 



Fig. 4: Observed major axis position-velocity diagram of NGC 973 (black and gray contours) overlaid with different 
models. The green contours show the final (B) model. Other models are shown as red contours. Contour levels are -2.3, 
2.3 (1.5 (t), 7.0 and 9.3 mJy beam“^. The greyscale corresponds to the observations. 


small warps of only a few degrees (as would be the case 
here), these changes are negligible. 

Fig. |A.l| of the appendix compares the model data cube 
to the observations. The satellite galaxy is indicated by ar¬ 
rows in the relevant channel maps. The main parameters 
describing the final model are given in Table and Fig. 
[29] The inclination and the position angle have errors of 
2.0° and 3.0° respectively. The scale height has a value of 
0.2lg'2 kpc. We note that this result seems to be consistent 
with a scale height of 0 kpc, which is of course unphysical. 
Unfortunately, due to the large observational beam it is not 
possible to determine a meaningful lower limit for the scale 
height of NGC 973. 

4.2. UGC4277 

With a distance of 76.5 Mpc, UGC4277 is the most dis¬ 
tant of the FIEROES galaxies. It is classified as an Sc or 
Scd galaxy and, like NGC 973, the details of its atomic gas 
content have not been investigated before. 


4.2.1. The data 


As can be seen in Eig.|^ both the stellar disk and the atomic 
gas disk of UGC4277 are slightly warped. Similarly to the 
case of NGC 973, the limited observing time results in a 
rather low signal-to-noise ratio in our data cube, making 
the individual channel maps quite noisy and the shape of 
the emission quite irregular (see Fig. A.2 of the appendix). 
This is also reflected in the major axis position-velocity di¬ 
agram (see Fig. later), where chunks of emission seem 
to be missing. However, the higher brightness contours, 
which should be less affected by the noise, are also very 
irregular. This suggests that UGC 4277 has a rather clumpy 
Hi disk. Similarly, the FIR morphology, especially in the 
PACS bands, suggests that the dust disk of UGC 4277 is 
also rather irregular ( jVerstappen et al.|2013 1. 

From our zeroth moment map we measure a total Hi 
flux of 14.7 Jy km s“^, corr esponding to an atomic hydro¬ 


gen mass of 2.03 x 10^° Mq. Huchtmeier & Richter ( 
report a single dish flux of 1 7.9 Jy km s~^, which was cor- 


( |1989| 


rected to 18.9 Jy km s ^ by Springob et al. (2005 1 . 



RA (J2000) 


Fig. 5: Top: False colour TNG V-band image of UGC 4277 
from [Bianchil ( |2007] l. Bottom: Hi contours overlaid on the 
same V-band image. Contours start at 3.61 x 10^° atoms 
cm“^ and increase as 2, 4, 6, 8 and 10 times this value. 
North is up. East is to the left. 


4.2.2. Models 

Eollowing the standard procedure we first constructed a ba¬ 
sic model consisting of a single disk with a constant incli¬ 
nation and scale height. This B model has an inclination 
of 87.3° and a rather large scale height of 1.8 kpc (but 
see last paragraph of this section). In the slightly warped 
outer edges of the disk the position angle deviates at most 
3.6° from the value in the inner disk. This model matches 
the observed total Hi map rather well (Eig. [^, although it 
does not reproduce the highest brightness contours. If we 
compare the model channel maps (Eig. [7|) and major axis 
XV-diagram (Eig. to the observations, we see the same 
problems as we encountered for NGC 973. The model often 
does not match the exact shape of the observed emission 
in the channel maps and lacks the observed high surface 
brightness contours at intermediate velocities. 
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Fig. 7: Representative channel maps from the observed data cube of UGC4277 and the various models. Contour levels 
are -1.5, 1.5 (1.5 (t) and 4.0 mJy beam“^. The black contours show the observations, with negative contours as dashed 
grey. The green contours represent the final (B) model. Other models are shown as red contours. The systemic velocity is 
5370±6kms-i. 


We therefore followed the same modelling strategy and 
constructed a B2 model by imposing a constant inclina¬ 
tion of 90°. TiRiFiC found the optimal corresponding scale 
height to be 2.2 kpc. As can be seen in Fig. this model 
does not solve the problems mentioned in the previous para¬ 
graph. The higher contours at intermediate velocities are 
still not reproduced, and additionally the contours in the 
central region of the XV-diagram are now overestimated. 
In the channel maps in Fig.|^this model also brings no im¬ 
provement. Finally we note that the zeroth moment map 
of the B2 model does now include the highest contours 
from the observed moment-0 map (Fig. |^, but they are lo¬ 


cated too far towards the center of the galaxy. This is actu¬ 
ally the direct consequence of underestimating the emission 
at intermediate velocities and overestimating the emission 
close to the systemic velocity, as we described for the XV- 
diagram. 

Next, we again manually scaled up the surface bright¬ 
ness distribution of the B2 model to match the high con¬ 
tours at intermediate velocities. In the end an increase of 
25 percent on the approaching side and 35 percent on the 
receding side was needed to more or less achieve this {B3 
model). However, from Fig. [^it is clear that the shape of 
the model contours is quite different from the observed 




































B: Basic 


B2: IncI = 90' 


B3: Extra gas at rec. side 



Fig. 8: Observed major axis position-velocity diagram of UGC 4277 (black and gray contours) overlaid with different 
models. The green contours represent the final (B) model. Other models are shown as red contours. Contour levels are 
-1.5, 1.5 (1.5 (t), 3.0 and 6.0 mJy beam“^. The greyscale corresponds to the observations. 
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of the irregular jumps of the latter), there is no reason to 
alter the geometry of the model and we take the B model 
as our final model. The observed data cube is compared to 
the model cube in Fig. |A.2| of the appendix and the main 
parameters of the model are shown in Fig. and listed in 
Table The values of the inclination, position angle and 
scale height have an uncertainty of 1.5°, 2.0° and 1.0 kpc 
respectively. 

The scale height of the final model is remarkably high. 
However, the large uncertainty on this value indicates that it 
should not be trusted blindly and still allows a lower, more 
common value of the scale height as well. In addition, a flar¬ 
ing gas disk, with a large scale height in the outer regions 
but a significantly lower scale height in the inner disk, is 
also not in disagreement with the data. The low signal-to- 
noise ratio and irregular shape of the observed emission in 
the channel maps unfortunately make it impossible to dis¬ 
tinguish between these options and put more stringent con¬ 
straints on the scale height. 

4.3. 1C 2531 


Fig. 6: Total Hi maps of the various models discussed here 
as compared to the observed total Hi map of UGC 4277. 
The B model is our best model. Contour levels are the same 
as in Fig.|^ 


contours, which are clumpy and follow an irregular pat¬ 
tern. Additionally, the model now strongly overestimates 
the emission near the systemic velocity and at the terminal 
edges. We note that the innermost contour of the B3 model 
between 5540 and 5620 km s“^ is actually a hole in the 
3CT-contour and not a 6cr-contour. The agreement with the 
observations in the total Hi map (Fig. and in the chan¬ 
nel maps (Fig. |7]l is also clearly a lot worse than for the B 
model. 

Similarly to NGC 973, the combination of a low signal- 
to-noise ratio and an irregular gas disk makes that little can 
be done to improve the model. Since, again, the shape of the 
model emission in the channel maps does not differ from 
that of the observed emission in a systematic way (on top 


IC2531 is the only Southern galaxy in our sample and is 
located at a distance of 36.8 Mpc. It is classified as an Sb or 
Sc type galaxy and has a prominent peanut-shaped bulge. 
Patsis & Xilouris (2006 1 ascribe this to a barred potential, 
although optical emission line observations do not show 
the characteristic line splitting in the bulge region ([Bureau 


& Freeman||1997[|1999[ Kregel & van der Kruit|2()04| l. A 

false colour V-band image taken with the FauUces Telescope 
South is shown in the top panel of Fig. The stellar disk 
shows a slight warp and is more extended on the approach¬ 
ing (NE) side, where it is also slightly misaligned with the 
dust lane. On the receding (SW) side, the dust lane shows 
a weak bump. The Hi content of 1C 2531 was already stud¬ 
ied in the series of papers by Kregel et al. ( Kregel et al. 
2004 Kregel & van der Kruit|20()4 1. They focused mainly 
on deriving the Hi rotation curves, with the ultimate goal 
to analyse the contribution of the stellar disk to the (inner) 
rotation curve of a spiral galaxy. IC2531 was also part of 
the sample of O’Brien et al. ( 2010d|b|a|c| l and later Peters 
et al. ( 2013| l, who studied the rotation and flaring of the Hi 
disks of a number of edge-on galaxies to derive the three- 
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Fig. 9: Top; False colour V-band image of 1C2531 taken 
with the Faulkes Telescope South. The data were cleaned 
and combined using standard data reduction techniques. 
Bottom: Hi contours overlaid on the same V-band image. 
Contours start at 2.98 x 10^° atoms cm“^ and increase as 
4, 8, 12, 16, 20 and 24 times this value. North is up. East is 
to the left. 



dimensional shape of their dark matter halos. These stud¬ 
ies, however, did not yet include the possibility of a ra¬ 
dially varying inclination. Kregel et al. ( |2004| l discovered 
two companion galaxies at projected distances of 12' (128 
kpc) and 14' (150 kpc) and with velocity separations of 200 
km s“^ and 182 km s“^ in their study of the Hi content of 
1C 2531. 
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Fig. 10: Total Hi maps of the various models discussed here 
as compared to the observed total Hi map of IC 2531. The 
W model is our best model. Contour levels are the same as 
in Fig. 1^ 


4.3.1. The data 

A contour plot of our total Hi map is shown in the bot¬ 
tom panel of Fig. overlaid on the V-band image from the 
top panel. Like the stellar disk, the Hi disk is also slightly 
warped and more extended towards the NE side. This was 
already observed by [Kregel et al.| ( |2004| l. The individual 
channel maps show that the behaviour of the Hi disk some¬ 
what deviates from that of a perfect regularly rotating disk 
(see Eig. |A.3| in the appendix). The approaching and reced¬ 
ing side behave quite differently, for example regarding the 
shape of the channel maps and the (projected) peak inten¬ 
sities. In addition, (projected) extraplanar features seem to 
appear and disappear irregularly throughout the cube. 


We measure a total Hi flux of 42.7 Jy km 


c-l 


which 

corresponds to a total Hi mass of 1.37 xlO^'^ Mq at the 
assumed distance of 36.8 Mpc. Our total Hi flux is remark¬ 
ably higher than the flux reported by O’Brien et al.| ( |2010d| l 
(13.1 Jy km s~^), but in relative ly good agreement with the 
values from 


Kregel et al.| ( [2004| (44. IJy km s" ^) and|P etersj 


et al. (20131 (41.7 Jy km s ^). These studies are based on 


the same data set as the current analysis with the exception 
of [Kregel et al. ( [2004[ l, who did not include the observations 
from January 2002 taken in the 750A configuration (see our 
Table ^and their Table 1). The single dish 21-cm fluxes re¬ 
ported on NED range fro m 33.1 Jy km s~^ (Koribalski et al. 
|2004[ ) to 52.3 Jy km s“^ ( Mathewson & Eord|1996 1 . 


4.3.2. Models 

Imposing the basic geometry, the best agreement of the B 
model with the data was achieved with an inclination of 
87.5° and a scale height of 1.1 kpc. The warp deviates at 
most 4° from the inner disk. Comparing the total Hi map 
of the basic model of IC2531 to the observed map (Eig. 
10 1 , it is clear that the low-density contours of the model 
disk are too thin. An obvious explanation for this would be 
that the scale height is underestimated or the inclination is 
overestimated, or a combination of both. In that case, the 
model emission (or at least a part of it) in the individual 
channel maps would be systematically thinner than the ob¬ 
served emission. Eig. [TT[ shows that this is not the case. The 
apparent thickness of the disk in the observed moment-0 
map is in fact caused by the extraplanar features indicated 
in the previous section. At the sensitivity of our data these 
do not seem to be part of a regular structure, and it is not 
possible to model them in the frame of an axisymmetric 
disk. We therefore excluded these features in the modelling 
and focused on the underlying disk. As a consequence, all 
models appear slightly too thin in the total Hi map. Deeper 
observations might reveal the true nature of these irregular 
features and whether they are part of an extended, low col¬ 
umn density structure such as a thick disk or a large-scale 
line of sight warp. 
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Fig. 11: Representative channel maps from the observed data cube of IC2531 and the different models. Contour levels 
are -1.6, 1.6 (l.Scr) and 6.5 mJy beam“^. The black contours show the observations, with negative contours as dashed 
grey. The green contours represent the final (W) model. Other models are shown as red contours. The dashed blue lines 
illustrate the triangular-like shape of the observed emission. The systemic velocity is 2455 ± 6 km s“^. 


Ignoring the extraplanar emission, the comparison of 
some key channel maps in Fig. [TT] particularly on the re¬ 
ceding side, shows that the geometry of the basic model is 
wrong. Indeed, in the observed channel maps we see that 
the emission has a triangular-like shape, with the inner part 
(i.e. the part closest to the center of the galaxy) being a bit 
thinner than the outer part. In addition, the outer tip of the 
observed emission generally appears quite flat, i.e. it looks 
more like a straight vertical line than a round arc. This is 
illustrated in Fig. by the blue dashed lines in the map 
at 2530 km s“^, where the effect is most obvious. The B 
model, on the other hand, has a round shape. The differ¬ 


ence between both might be difficult to see in Fig. but 
is obvious in Fig. The triangular shape of the observed 
emission indicates the presence of a line of sight warp or 
a flare. We therefore continued by investigating these two 
geometries. 

The best agreement with the inner parts of the observed 
emission in the channel maps was obtained with an incli¬ 
nation of 89.5° and a scale height of 1.0 kpc for the inner 
disk. A flaring {F) model was constructed by gradually in¬ 
creasing the scale height in the outer rings while fixing the 
inclination to 89.5°. The flare starts at a radius of 25 kpc on 
both sides of the disk and grows to a maximum scale height 
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Fig. 12: Channel map at 2530 km s ^ from the observed cube of IC 2531 (black and grey contours), the basic (B) model 
(red contours) and the final (W) model (green contours). Contour levels are -1.6, 1.6 (1.5cr) and 6.5 mJy beam“^. The B 
model is deary too round. The W model has a more triangular shape and provides a better fit to the data. 


of 1.3 kpc on the approaching side and 1.8 kpc on the re¬ 
ceding side. Fig. [TT] shows that the flaring model provides 
a better fit to the data on the inner side of the emission, but 
the outer edge of the model is still systematically too round 
on the receding side. This rounding is inherent to a flaring 
geometry and can not be removed by changing the details 
of the flare. 

A model with a warp along the line of sight naturally 
explains these issues. The best fitting model (W) has a con¬ 
stant scale height of 1.0 kpc, a central inclination of 89.5° 
and a line of sight warp that reaches an inclination of 85° on 
the receding side and 87° on the approaching side. Figures 
11 and 12 show that the triangular shape of the emission 


on the receding side is now nicely matched by the model. 
The high-density contours are a bit underestimated in some 


channels (e.g. at 2569 and 2609 km s ^ in Fig. Ill, but 


4.3.3. Spiral arms 

The terminal side of the major axis position-velocity dia¬ 
gram (Fig. [T3 ] i is nicely reproduced by the W model. This 
indicates that we recover the rotation curve and the global 
Hi distribution of the disk quite well. Not reproduced by our 
model is the high-density ridge that is clearly visible along 
the low-velocity edge of the diagram. This ridge is present 
on both sides of the galaxy, albeit at lower contour levels on 
the approaching side, and was already noticed by |Kregel| 


et al.| (2004 1 on the receding side. After inspecting their 


this is most likely due to the presence of a spiral arm and 
will be discussed below. The agreement with the data on 
the approaching side seems less good in Fig. 0 although 
it is still better than for the previous models. However, the 
misalignment between the stellar disk and the dust lane, the 
larger extent seen in the optical and in HI, and the irregu¬ 
lar channel maps all suggest that this side of the galaxy is 
somewhat disrupted and not simply a smooth disk. The aim 
of this modelling is therefore not to obtain an exact match 
in every channel map, but to constrain the main properties 
of the Hi disk. 

Rather than just one specific type of geometry, real 
galaxy disks usually contain a combination of different fea¬ 
tures. Therefore we also investigated the combination of a 
line of sight warp with a flare. A substantial line of sight 
warp is in any case required to reproduce the flattening of 
the outer edge of the emission in the channel maps. Adding 
to this a flare of the same order as in the F model makes 
the emission too thick in the outer channels. Adding a more 
modest flare does not cause significant changes compared 
to the W model. This is shown for a flare with a maximum 
scale height of 1.3 kpc on both sides as the WF model in 
Fig. 11 Since the addition of a flare does not result in a 


siginificant improvement of the model, but does introduce 
extra free parameters, we discard this additional feature and 
stick to the simple line of sight warp geometry {W model). 


full-resolution data cube, they conclude that the ridge is a 
superposition of multiple rings. On the other hand, in the 
next paper of this series Kregel & van der Kruit ( 2004| l ar¬ 
gue that the coincidence of HII emission with the Hi ridges 
in many galaxies could indicate the presence of spiral arms. 

We tested the first hypothesis by modelling the ridge as 
a set of rings in different ways. A first possibility is that 
the ridge is in fact a projection along the major axis of 
some high-density rings at large radii, rotating at the nor¬ 
mal terminal velocities. In practice this approach requires 
increasing the densities of the outer rings. On the approach¬ 
ing side, we were unable to construct a ridge in this way, 
no matter how much the densities were increased. On the 
receding side, this strategy led to a huge overestimation 
of the intensities in the tip of the XV-diagram for the W 
model. For a flaring model with a constant inclination of 
88°, the outer rings are more strongly projected along the 
major axis. Hence a smaller increase of the outer densities 
suffices to create a ridge in the XV-diagram. But even in 
this case the intensities in the tip of the XV-diagram were 
overestimated. 

A second, more artificial possibility would be that the 
ridge is a superposition of multiple rings at different radii, 
rotating at velocities considerably lower than the terminal 
velocities. Such a model, based on the W model, is shown 
in Fig.[^ Although it is indeed possible to reconstruct the 
observed ridge, this approach requires adding a large num¬ 
ber of low-velocity rings on top of the main disk (see Fig. 
[T4|for the details of these rings) and is highly unphysical. 

Based on our modelling we thus exclude the interpreta¬ 
tion of the high-density ridge as a superposition of rings. A 
prominent spiral arm provides a much more natural expla¬ 
nation. Unfortunately, the edge-on orientation offers little 
constraints and modelling the spiral arms would be highly 
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Fig. 13: Observed major axis position-velocity diagram of IC2531 (black and gray contours) overlaid with different 
models. The green contours represent the final (W) model. Other models are shown as red contours. Contour levels are 
-1.6, 1.6 (l.Scr), 4.9 and 8.1 mJy beam“^ and the greyscale corresponds to the observations. The high-density ridge along 
the low-velocity edge on both sides (indicated by the arrows) is not reproduced by the W model. 



Fig. 14: Surface densities (top) and rotational velocities 
(bottom) of the rings added to the W model of IC2531 in 
order to reproduce the high-density ridge in the major axis 
XV-diagram. Blue indicates the approaching side, red the 
receding side. 


speculative at best. We therefore do not attempt to include 
spiral arms in our model. Because of this, the peak densities 
in the total Hi map and in some channel maps are slightly 
underestimated by the W model. To prove that the latter 
is not in fact caused by a structural mistake in the model, 
we constructed a toy model {W+A) that reproduces the Hi 
ridge in the XV-diagram by locally enhancing the surface 
brightness in several arcs within the disk (Fig. [D). Figures 
ITOl and fTTI show that with this addition the model indeed 
reproduces the observed peak densities. 

Despite the improved agreement with the data at high 
column densities, adding arcs to a model without a clear 
physical motivation is a dangerous practice and not unam¬ 
biguous. The W+A model therefore remains a toy model 
and we take the W model as our final model for the atomic 
gas disk of IC 2531. The total Hi map of the W model is 
compared to the observed map in Fig. 10 Our final model 
data cube is compared to the observed cube in Fig. A.3 of 
the appendix. The main parameters of the final model are 
given in Fig. 29 and Table The errors on the inclination, 
position angle and scale height are 1.3°, 1.8° and 0.5 kpc 
respectively. 


4.4. NGC4217 

NGC4217 is located at a distance of 19.6 Mpc and is classi¬ 
fied as an S Ab galaxy. Its atomic gas content has previously 
been studied by [Verheijen & Sancisi| ( |2001| l; |Verheijen| 
( |2001) l, who used the kinematics of the atomic gas in a sam¬ 
ple of galaxies in the Ursa Major cluster to investigate the 
statistical properties of the Tully-Fisher relations. However, 
their sample consisted largely of moderately inclined galax¬ 
ies and their method to derive the surface density distribu¬ 
tion of the Hi is actually not suitable for an edge-on geome¬ 
try, with Verheijen and Sancisi stating that the “method for 
extracting the surface density profiles from integrated Hi 
maps breaks down for nearly edge-on systems.” 


4.4.1. The data 


In the total Hi map of NGC4217 (Figs. [H® two blobs 
of gas on either side of the main disk immediately catch the 
eye. They suggest the presence of a ring, which, strikingly, 
has no optical counterpart in the r-band image in Fig. 

In what follows we will refer to this feature as the outer 
ring. Some filaments can also be seen extending out of the 
plane of the galaxy. The Hi disk is not significantly warped, 
apart from a slight upward bend (i.e. towards the NW in Fig. 
0 where it meets the outer ring. This upward bend is also 
visible in the stellar and dust disks. 


If we separate the outer ring from the main disk in our 
data (see section |4.4.2| for a discussion on the nature of the 
outer ring), we measure a total Hi flux of 27.6 Jy km s“^ 
for the main disk and 4.2 Jy km s“^ for the outer ring. At 
a distance of 19.6 Mpc, this translates to atomic hydrogen 
masses of 2.50x10® Mq and 3.81x10® M©. It should be 
noted that in the central channels of the data cube the emis¬ 
sion from the outer ring and the main disk are projected 
on top of each other, making it very difficult to separate 
both components. Our estimate of the mass of the outer ring 
only includes the emission that unambiguously belongs to 
this component, so it should be regarded as a lower limit. 


Verheijen & Sancisi (2001 1 also analysed the Hi content of 


NGC4217, based on the same data as used here, but did 
not make a distinction between the main disk and the outer 
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Fig. 15: Top: False color r-band image of NGC4217 (SDSS 


J121551.51 +470538.1) from the EFIGI catalogue ( jBaillard 
|et al.|20lT| l, based on data from the SDSS DR4. Bottom: Hi 
contours overlaid on the same r-band image. Contours start 
at 3.76 X 10^° atoms cm“^ and increase as 2.5, 4, 8, 12 and 
16 times this value. North is up, East is to the left. 


ring. They found a total Hi flux of 33.8 Jy km s“^, which 
is slightly higher than our total flux of 31.8 Jy km s“^. The 
reason for this is probably that we were more strict in sep¬ 
arating the real emission from the noise. (Corrected) single 


dish fluxes of 25.8 Jy km s ^ and 31.0 Jy km s ^ are re- 


-1 


ported by Springob et al.| ( 2005] l and Huchtmeier & Richter 
(19891 respectively. 


4.4.2. Models 


In the major axis position-velocity diagram (see Eig. 18 
later) the outer ring corresponds to the sudden drop in rota¬ 
tional velocities at major axis offsets of about ± 3 arcmin. 
The outer ring thus seems both spatially and kinematically 
separated from the main disk, indicating that it probably 
has an external origin. Because of this, we start by mod¬ 
elling only the main gas disk and add the outer ring as an 
extra feature in a later stage. 

An initial fit (B) imposing the basic geometry resulted 
in an inclination of 86.2° and a scale height of 0.5 kpc. Eig. 
17 compares this basic model to the observations. Ignoring 
the outer ring, the model already provides a good fit to the 
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Fig. 16: Total Hi maps of the various models discussed here 
as compared to the observed total Hi map of NGC4217. 
The B model is our best model. Contour levels are the same 


as in Fig. 15 


data. The only exception to this are four channels on the 
receding side, between 1059 and 1109 km s“^, and (to a 
lesser extent) two channels on the approaching side, be¬ 
tween 928 and 944 km s“^, where the model is significantly 
thicker than the observed emission in the inner regions. Two 
of these channels are shown in Fig. 17 at 928 and 1076 
km s“^, where the discrepant regions are indicated with ar¬ 
rows. This discrepancy could indicate a flaw in the model, 
but the fact that it only shows up in a few channels and that 
the affected channels are different on both sides (the cen¬ 
tral channel is at 1026 km s“^) indicates that it is probably 
just caused by local irregularities in the disk and does not 
reflect a large-scale geometric feature. Nonetheless, we in¬ 
vestigate whether a more complex geometry could provide 
a better fit. 

With the parameters of the basic model as inputs we 
conducted a new fit with TiRiFiC. The inclination and scale 
height were now allowed to vary in blocks of five rings in 
order to reveal a potential flare and/or line of sight warp, 
but at the same time smooth out the unphysical jumps that 
TiRiFiC usually finds. Fits in blocks of three and seven 
rings were also performed to verify the results. The incli¬ 
nation on the receding side was found to roughly increase 
towards the center, reaching a maximum of 89.5° in the 
central rings. For the inclination on the approaching side 
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Fig. 17: Representative channel maps from the observed data cube of NGC4217 and the different models of the main 
disk. Contour levels are -1.5, 1.5 (1.5tT) and 8.0 mJy beam“^. The black contours show the observations, with negative 
contours as dashed grey. The green contours represent the final (B) model. Other models are shown as red contours. The 
systemic velocity is 1022 ± 6 km s“^. The arrows highlight regions where the B model is significantly thicker than the 
observed emission, as explained in the text. 


and the scale height on both sides, the radial profile found 
by TiRiFiC was still dominated by unphysical jumps and 
did not display a clear trend. In the absence of strong con¬ 
straints we constructed a model with a line of sight warp by 
applying the rising trend of the inclination that was found 
for the receding side to the approaching side as well and 
imposing a constant global scale height. A new TiRiFiC fit 
found a value of 0.6 kpc for the latter. As can be seen in 


Fig. 17 the difference between this model (W) and the B 
model is very small and the discrepancy mentioned in the 
previous paragraph is still present. Moreover, looking at the 


r-band image in Fig. 15 a central inclination of 89.5° seems 
too high. Indeed, the bulge is clearly a lot stronger on the 
NW side of the galaxy (‘above’ the major axis) while the 
dust lane is located more towards the SE, suggesting that 
the galaxy is not seen perfectly edge-on. For these reasons 
we reject the W model and continue with a constant incli¬ 
nation. 

In a final attempt to model the thinning of the disk we 
added a flare to the B model. The F model shown in Fig. 


17 has a constant scale height of 0.1 kpc out to a radius of 


2.6 kpc, which then linearly increases to reach a maximum 
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Fig. 18: Observed major axis position-velocity diagram of NGC4217 (black and gray contours) overlaid with different 
models. The green contours show the final (B) model. Other models are shown as red contours. Contour levels are -1.5, 1.5 
(1.5 ct), 4.5, 9.0 and 13.5 mJy beam“^. The greyscale corresponds to the observations. We interpret the features indicated 
by the arrows as the projection of a clumpy ring along the major axis. 


value of 0.6 kpc in the outer regions of the disk. The dis¬ 
crepancy in the six channels mentioned earlier is still not 
completely resolved, although it is (slightly) less severe. 
However, at the same level as this flaring model improves 
the fit in these six channels, it is now too thin in many other 
channels. Adding a flare to the model therefore introduces 
many free and poorly constrained parameters without sig¬ 
nificantly improving the quality of the fit. We thus decided 
not to include a flare and to stick with the basic geometry (B 
model). It should be noted, though, that the difference be¬ 
tween the various models we discussed here is very small 
and hence the geometry of the gas disk is not very well 
constrained by our modelling. However, this does not sig¬ 
nificantly affect the other parameters. Between the differ¬ 
ent geometries investigated here, only the surface bright¬ 
ness distribution needed slight adjustments and the changes 
were of the order of 10% at most. 


4.4.3. The outer ring 


The next step in the modelling process is to determine the 
nature of the outer ring. We start by modelling it as a simple 
extension of the main disk, i.e. using the same inclination 
and scale height as in the disk. This model (model B+Rl in 
Fig. [19)1 performs well in the very outer channels, but in the 
other channels the projection of the outer ring along the ma¬ 
jor axis is clearly in disagreement with the data, especially 
on the approaching side. Arrows highlight this effect in the 
channel maps at 895, 944 and 1158 km s“^ in Fig. 19 This 
is also visible in the XV-diagram (Fig. 18 i and indicates that 
either the outer ring is not a complete ring, or that the in¬ 
clination should be lower, placing less gas along the major 
axis. We test the latter option in the B+R2 model. Fig. 
shows that this model does indeed solve the aforementioned 
projection problem, but cannot account for the clumpy fea¬ 
tures along the low-velocity edges of the XV-diagram (in¬ 
dicated with arrows for the B model). Moreover, the lower 
inclination makes this model too thick, as can be seen in the 
moment-0 map (Fig. [Tbl l and in the central channels in Fig. 
19 Finally, both these models also require a highly unphys¬ 
ical rotation curve on the receding side, with a rather abrupt 
decrease of about 60 km s“^ (see Fig. 291. All these argu¬ 


ments indicate that the outer ring is not simply an extension 
of the main disk. A much more natural interpretation would 
be that we are observing the accretion of a roughly coplanar 
satellite galaxy that has been gradually ripped apart, form¬ 
ing a clumpy arc or ring-like structure around NGC4217. 
This would explain the spatial and kinematical offset with 
respect to the main disk and the absence of a smooth pro¬ 
jection of the ring along the major axis. 


Modelling the outer ring as a set of individual clumps is, 
in principle, possible with TiRiFiC, but this exercise would 
not contribute much to the overall quality of the model 
and is beyond the scope of this paper. We therefore take 
the B model as our final model for the atomic gas disk of 
NGC4217 and note that a coplanar ring-like structure just 
outside of the disk is not included in this model. A compar¬ 
ison between the final model and the observed data cube is 
shown in Fig. A.4 of the appendix. The major axis position- 
velocity diagrams are compared in Fig. 18 and the main pa¬ 
rameters of the final model are shown in Fig. and Table 
The inclination and position angle have respective uncer¬ 
tainties of 1.5° and 2.0° and the error on the scale height is 
0.2 kpc. 


4.5. NGC5529 


NGC 5529 is an Sc type galaxy located at a distance of 49.5 
Mpc. Like IC2531, NGC 5529 was also part of the sample 
of galaxies analysed by |Kregel et al. ( |2004[ ); Kregel & van 
der Kruit (|2004|l. In the region of the peanut-shaped bulge. 


Kregel & van der Kruit (2004 1 find a typical figure-of-eight 


pattern in their optical emission line position-velocity dia¬ 
gram, indicating the presence of a bar. In the outer regions 
of the galaxy, the optical emission line profiles are again 
double peaked. NGC 5529 is part of a rich galaxy group 
with at least 16 other members ( |Irwin et al.|2007[ ). Its stel¬ 
lar disk shows a clear warp and extends out to larger (pro¬ 
jected) radii on the NW side. In their study of the atomic 
gas in this galaxy, Kregel et al. ( 2004| l discovered that 
NGC 5529 is connected to two of its companions, MGC 
H-06-31-085a and NGC 5529 B, through Hi bridges. 
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Fig. 19: Representative channel maps from the observed data cube of NGC4217 and the different models that include 
the outer ring. Contour levels are -1.5, 1.5 (1.5cr) and 8.0 mJy beam“^. The black contours show the observations, with 
negative contours as dashed grey. The green contours represent the hnal (B) model. Other models are shown as red 
contours. Because we want to highlight different features, the channels shown here are different from those in Fig.[^ The 
arrows highlight the projection along the major axis of the outer ring in the B+Rl model, which is clearly in disagreement 
with the data. 


4.5.1. The data 


Derived from the same WSRT data as used by Kregel et al. 
( [2004[ ), our total Hi map (Fig. 201 is very similar to theirs. 
On the approaching (NW) side the atomic gas follows the 
same moderate warp as the stars. On the receding (SE) side 
the Hi appears more strongly warped and shows a strong 
upturn beyond the stellar disk. The channel maps (see Fig. 


A.5 in the appendix) reveal that the Hi disk of NGC 5529 


is surrounded by a significant volume of extraplanar gas. 


especially on the approaching side, as a consequence of the 
interaction with the satellite galaxies. This will be discussed 
in further detail in the next section. 

It is not possible to unambiguously separate the actual 
Hi disk of NGC 5529 from the extraplanar gas and the satel¬ 
lite galaxies in the channel maps. In order to determine 
the total Hi mass of the galaxy, we took a conservative 
approach and neglected only the emission that unambigu¬ 
ously comes from the satellite galaxies. From the resulting 
moment-0 map we measure a total Hi flux of 46.5 Jy km 
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Fig. 20: Top; False colour V-band image of NGC5529 
taken with the 1.3m telescope at the Skinakas Observatory 
in Crete. The data were cleaned and combined using stan¬ 
dard data reduction techniques.. Bottom: Hi contours over¬ 
laid on the same V-band image. Contours start at 9.94 x 
10^® atoms cm“^ and increase as 4, 10, 20, 30, 40, 50 and 
60 times this value. North is up. East is to the left. 


s“^. This is slightly higher than the value of 43.2 Jy km 
s“^ found by Kregel et al. (2004 1 . The small difference be¬ 
tween both fluxes is most probably due to the rather subjec¬ 
tive distinction between the emission from NGC5529 and 
that from the satellites. At a distance of 49.5 Mpc our Hi 
flux corresponds to a total Hi mass of 2.69 x 10^° M©. 
For NGC 5529, [Huchtmeier & Richter] ( |1989| l report total 
HI fluxes of 26.8, 37.1 and 40.8 Jy km s obtained with 
different single dish telescopes. 
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Fig. 23: Total Hi maps of the various models discussed here 
as compared to the observed total Hi map of NGC 5529. 
The W2R model is our best model. Contour levels are the 


same as in Fig. 20 


4.5.2. Models 


Imposing a simple, one-component disk to start the mod¬ 
elling, the best agreement with the data was achieved with 
an inclination of 86.4° and a scale height of 0.4 kpc. The ap¬ 
proaching side of the resulting B model is dominated by a 
gradual warp that starts early on in the disk, at a radius of 15 
kpc, and reaches a maximum deviation of 10° with respect 
to the central position angle. The receding side shows a 
more abrupt warp, starting at a radius of 37 kpc and rapidly 
increasing to a maximum deviation of 16°. Comparing the 
B model to the data in the channel maps (Figs. pTjp^ and 
the moment-0 map (Fig. 23 i, we see that the model already 
provides a reasonable approximation to the data, although 
there are still some issues to address. Firstly, part of the ob¬ 
served channel maps on the approaching side contain addi¬ 
tional emission that is trailing behind the model. This emis¬ 
sion is indicated with arrows at 2615, 2697 and 2780 km 


c-l 


in Fig. 21 and is also visible in Fig. A.5 between 2598 


and 2763 km s ^. However, in the outermost channels the 


model reproduces all the observed emission very well. In 
addition, the extra emission never exceeds the outer edge 
of the (integrated) model disk. These two arguments indi¬ 
cate that the model disk extends out far enough radially and 
is rotating at the correct velocity and that the extra emis¬ 
sion is not the projection of larger rings. We thus interpret 
this emission as extraplanar gas resulting from the interac¬ 
tion with MGC -1-06-31-085a and do not include it in our 
models of the main disk of NGC 5529. 

A second, more subtle deficiency of the B model is best 
visible in the innermost channel maps, for example at 2780 
and 2879 km s“^ in Fig.and is highlighted for the 2829 
km s“^ channel in Fig. |22| While the model emission in 
these channels is oriented more or less perpendicular to the 
major axis (vertical in Figs. [2T][22) i, the observed emission 
shows a slant and is oriented more diagonally, from the top 
right to the bottom left. This discrepancy cannot be resolved 
by changing the geometry of the disk. Instead, we model the 
slant by introducing inward radial motions of -15 km s“^ in 
the entire disk (BR model). We note that from our Hi data 
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Fig. 21: Representative channel maps from the observed data cube of NGC5529 and the various models. Contour levels 
are -0.8, 0.8 (l.Scr) and 5.6 mJy beam“^. The black contours show the observations, with negative contours as dashed 
grey. The green contours represent the final {W2R) model. Other models are shown as red contours. The systemic velocity 
is 2830 ± 6 km s“^. The arrows highlight extraplanar gas that probably results from the interaction with the companions. 


alone it is in principle not possible to determine which half 
of the disk (north or south of the major axis in Fig.[20|) is the 
near side and which is the far side. In the optical, however, 
we clearly see that the dust lane is mainly located north 
of the major axis, indicating that this is the near side. The 
radial motions in our model are inward under this assump¬ 
tion. The addition of a radial inflow significantly improves 
the agreement with the data, especially in the inner chan¬ 
nels. On the other hand, this purely kinematical adjustment 
has no influence on the moment-0 map. 


The total Hi map of the BR model is therefore identi¬ 
cal to that of the B model. In general it is in good agree¬ 
ment with the observed moment-0 map (Fig. [2^ , but at the 
lowest contour levels it is too thin. Adding a flare to the 
model can immediately be ruled out as potential solution 
to this problem. Indeed, in the outermost channel maps we 
see that the model emission is thick enough and matches 
the data very well. Adding a flare would make the model 
significantly too thick in these channels. A second alterna¬ 
tive consists of adding a line of sight warp. We tested this 
option by gradually lowering the inclination from 86.4° to 
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Fig. 22: Channel map at 2829 km s ^ from the observed data cube of NGC 5529 (black and grey contours), the B model 
and the BR model. Contour levels are -0.8, 0.8 (1.5tT) and 5.6 mJy beam“^. The orientation of the B model is clearly too 
vertical. Adding a radial inflow significantly improves the agreement with the data. 


82.0° between radii of 40 kpc (2.8') and 47 kpc (3.3') and 
keeping it fixed at 82.0° beyond that. The moment-0 map 
of the resulting WIR model now has the correct thickness, 
but in the channel maps this model does not provide a con¬ 
sistent improvement over the BR model (Fig.[2T]l. In some 
channels the observed emission broadens on the outer side 
and shows a triangular-like shape, which is nicely repro¬ 
duced by the WIR model (e.g. at 2615 km s“^ in Fig. 21 1 . 
In other channels, however, this broadening is less strong 
or even absent and the emission from the WIR model is 
too thick on the outer side (e.g. at 2879 and 2961 km s“^). 
The absence of a systematic trend in the data leads us to 
believe that the strong broadening of the observed emission 
in some of the channels is not a real feature of the disk, but 
is rather caused by the extraplanar gas from the interaction 
with the satellites. As a consequence, the lowest contours 
in the observed moment-0 map most probably also reflect 
the extraplanar gas and not the real thickness of the disk. 


Notwithstanding the arguments from the previous para¬ 
graph, the data do show indications of a small line of sight 
warp. Indeed, comparing the BR model to the data in the 
channel maps we see that the outer tip of the model emis¬ 
sion is generally a bit too extended and sharp. In order to 
correct for these effects without affecting the central chan¬ 
nels we constructed the W2R model. This has a central in¬ 
clination of 87.5° and a constant scale height of 1.1 kpc. 
The inclination then gradually decreases to 85.0° between 
40 kpc and 47 kpc from the center and stays constant be¬ 
yond that. This model provides a subtle but significant im¬ 
provement over the BR model. 


4.5.3. Spiral arms 


The terminal edges of the major axis XV-diagram (Fig. 
24 left panel) are reproduced rather well by the W2R 
model, both on the approaching and the receding half 
Similarly to the case of IC 2531, however, the XV-diagram 
of NGC 5529 shows an additional Fll ridge that is not re¬ 
produced by our model of the gas disk (indicated by ar¬ 
rows in the left panel of Fig. [24|. An d similarly to the case 
of IC2531, Kregel & van der Kruit ( |2004 1 found that this 
ridge coincides with a ridge of Ha emission, suggesting a 
spiral arm. Given the many uncertainties that are involved, 
we made no attempt to include spiral arms in our model, but 
we did again construct a toy model {W2R+A) that repro¬ 
duces the Hi ridge by locally enhancing the surface bright¬ 


ness in various arcs in the disk. The major axis XV-diagram 
of this model is shown in the right panel of Fig. and the 
corresponding zeroth moment map and channel maps can 
be seen in Figs. [^and[2^ We see that with the addition 
of the arcs, the peak intensities of the model are indeed in 
agreement with the data. 


As for IC 2531, we note that the addition of arcs to the 
W2R model is not really physical and not unambiguous. We 
therefore consider the W2R+A model as a toy model and 
take the W2R model as our final model for the atomic gas 
disk of NGC 5529. A comparison of this model to the full 
observed data cube is shown in Fig. |A.5| and the major axis 
position-velocity diagrams are compared in Fig. 24 Table 
[^and Fig.[29]show the main parameters of the final model. 
The inclination, position angle and scale height have errors 
of 2.0°, 3.0° and 0.7 kpc respectively. 


4.5.4. About the radial motions 


The central channel maps of the observed cube show a clear 
slant that indicates the presence of radial motions in the 
disk. To model this effect we used a global inflow in the en¬ 
tire gas disk. We acknowledge that such a parametrisation 
is most probably too simplistic, but the data do not provide 
sufficient constraints to model the radial motions in more 
detail. In order to reproduce the observed slant, an inflow 
of -15 ±5 km s“^ is definitely required in the outer half of 
the disk. On the other hand, radial motions in the inner parts 
of the disk only cause very subtle changes in the channel 
maps. Unfortunately these subtle changes are completely 
dominated by the effects of the extraplanar gas in the ob¬ 
served cube. It is therefore not possible to determine from 
the data whether an inflow is also present in the inner disk or 
not. In fact, even a small outflow in the inner regions is not 
inconsistent with the data. As an example of this. Fig. 
compares our final model to a similar model {W2R2) with 
no radial motions in the inner disk. The difference between 
both is indeed very small and it is not possible to determine 
which model is better. In order to minimise the number of 
free parameters we therefore opted to use a global inflow in 
the entire disk. 


The velocity of the radial motions derived in t his work 
is comparable to the inflow of -20 km s“^ found by de Blok 
et al. (2014|l in the outer part of the gas disk of NGC 4414. 


Zschaechner et al. (2012 1 report a strong, local inflow of ■ 
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Fig. 24: Observed major axis position-velocity diagram of NGC 5529 (black and gray contours) overlaid with different 
models. The green contours show the final W2R model, the red contours show the W2R+A model. Contour levels are -0.8, 
0.8 (1.5tT), 5.6, 11.2 and 16.8 mJy beam“^. The greyscale corresponds to the observations. The arrows indicate a Hi ridge 
that is not reproduced by the W2R model. 


50kms ^ in the inner part and an outflow of 10 km s ^ in 
the outer part of the gas disk of NGC 4565. 

How should we interpret this radial inflow? It was al¬ 
ready mentioned that NGC 5529 shows strong evidence for 
the presence of a bar and prominent spiral arms. These fea¬ 
tures are indeed known to generate non-circular flows in the 
disk of a galaxy. However, theory and simulations show that 
bars and spiral arms induce an inflow of gas inside the coro¬ 
tation radius and an outflow outside the corotation radius 
( [Kenney 1 1994[[Kormendy & Kennicutt|2004| l. Since coro¬ 
tation generally occurs at only several kpc from the center 
of a galaxy ( jElmegreen et al. 1992 1 , it is unlikely that the in¬ 
flow we And for NGC 5529 is driven by a bar or spiral arms. 
Instead it might be related to ongoing accretion of gas from 
the satellite galaxies. Indeed, when the accreted gas has a 
lower angular momentum than the gas in the disk, it can 
generate a radial inflow in the galaxy (Matteucci|2012l. 


4.6. NGC5907 


The closest member of our sample is NGC 5907, located at 
a distance of 16.3 Mpc. It is a disk-dominated galaxy with 
no significant bulge visible. In contrast to its morphologi¬ 
cal classification as an SA(s)c type galaxy, |Garcia-Burillo| 
et al. \\991 1 saw strong evidence for the presence of a bar 
in their CO(l-O) observations of the nucleus of NGC 5907. 
Almost 40 years ago, Sancisi ( [1976 1 already noticed that 
the Hi disk of NGC 5907 is strongly warped. This warp is 
also visible in the optical, although less obvious. NGC 5907 
is a member of the NGC 5866 group (or LGG 396), but the 
three other members of this group are located at large dis¬ 
tances. NGC 5907 was therefore long considered as the pro¬ 


totype of an isolated galaxy with a warped disk, until jShang 
et al. ( 1998|l discovered a small companion, PGC 54419, in 


their Hi map and an extended stellar tidal stream around the 
galaxy in their optical observations. Deeper optical imaging 
by Martlnez-Delgado et al. (2008 1 later revealed even more 


loops of stellar debris. The latter also conducted N-body 
simulations showing that most of the stellar loops they ob¬ 
served can be explained by the accretion of a single satellite 
galaxy. 


4.6.1. The data 

The total Hi map of NGC 5907 as derived from our data 
is compared to an r-band image from the SDSS in Fig.|25l 
In the optical image we see that the dust lane is located 
slightly below or west of the major axis. We therefore in¬ 
terpret that side as the front half of the disk and the east 
side as the back half. The Hi reveals a more complex struc¬ 
ture than would be expected from the optical image. Indeed, 
the main gas disk is strongly warped and slightly asym¬ 
metric, with the receding (NW) half showing a slight up¬ 
turn in the outer regions. Additionally, an extended body of 
low-density gas seems to surround the main disk (in projec¬ 
tion). Strikingly, this low-density gas is located almost en¬ 
tirely above (i.e., east of) the disk. Finally, the companion 
galaxy PGC 54419 is visible toward the NW of NGC 5907. 
We note that the atomic gas content of NGC 5907 might ac¬ 
tually extend even further than seen in our map since a sig¬ 
nificant amount of low-density gas is likely still hidden in 
the noise in our channel maps. Deeper WSRT observations 
with a total integration time of about 16 x 12 hours were 
recently conducted by G. Jozsa and a detailed modelling of 
these data is currently underway (Yim et al., in prep.). 

Excluding the companion, we measure a total Hi flux of 
317.6 Jy km s“^. This corresponds to an atomic hydrogen 
mass of 1.99 x 10^° Mq if we assume a distance of 16.3 
Mpc. NED lists a variety of single dish fluxes, r anging from 


36.6 Jy km s ^ ( 

Huchtmeier & Richter||l98S 

1 to 282.5 Jy 

km s ^ ( 

Springo 

j et al.|2005 i. Shang et al. ( 

19981 did not 


measure the total Hi flux of NGC 5907. Somewhat surpris¬ 
ingly, our interferometric flux measurement is higher than 
all the single dish fluxes. Keeping in mind the large angular 
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Fig. 25: Top: False color r-band image of NGC 5907 (SDSS 


J151553.77 +561943.5) from the EFIGI catalogue ( iBalUard 
|et al.|20lT| l, based on data from the SDSS DR4. Bottom: Hi 
contours overlaid on the same r-band image. Contours start 
at 1.01 X 10^° atoms cm“^ and increase as 5, 15, 30, 45, 60 
and 75 times this value. North is up. East is to the left. 


extent of NGC 5907, the limited HPB W of many of the sin¬ 
gle dish observations is likely part of the reason for this. In 
addition, uncertainties in the flux calibration generally also 
cause differences. 

4.6.2. Models 

The best fitting single-component (B) model has a constant 
inclination of 85.4° and scale height of 1.2 kpc. The strong 
warp on the approaching side is reproduced with a sudden 
strong decrease of the position angle at a radius of 26 kpc 
(see e.g. Eig.[26|). This decrease then becomes more gradual 
and continues out to the edge of the disk, where the position 
angle reaches a maximum deviation of 12.9° with respect 
to the center. On the receding side the position angle ini¬ 
tially shows a similar behaviour with a maximum deviation 
of 12.9° with respect to the center. However, at a radius of 
45 kpc the position angle then starts to increase again, even¬ 
tually deviating only 2.8° from the central value in the outer 
regions of the disk (Eig. p6| ). We also note that in the best 
fitting B model the position angles in the inner, flat part of 
the disk differ by 1.8° between the approaching and reced- 
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Fig. 26: Total Hi maps of the various models discussed here 
as compared to the observed total Hi map of NGC 5907. 
The WF3 model is our best model. Contour levels are the 


same as in Fig. 25 


ing sides. Comparing this model to the observed data cube 
we see that in the outer channels (e.g. at 424, 486, 548, 837 
and 898 km s“^ in Fig. 27 1 the model is clearly too thick in 
the central regions of the disk. Additionally, the inclination 
in the outer regions of the model still seems too high, since 
the observed emission in the channel maps clearly ‘opens 
up’ more than the model does. We address these issues one 
by one. 

Firstly we note that the central parts of the model disk 
appear too thick primarily in the outer channels, near the 
terminal velocities, where the inclination only has a neg¬ 
ligible effect. This indicates that the problem is related to 
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Data B: Basic F: Flare WF1: Warp + Flare 




Fig. 27: Representative channel maps from the observed data cube of NGC5907 and the various models. Contour levels 
are -0.5, 0.5 (1.5cr) and 6.4 mJy beam“^. The black contours show the observations, with negative contours as dashed 
grey. The red contours represent the models. The systemic velocity is 664 ± 7 km s“^. 
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Fig. 27: continued. The final (WF3) models is shown as green contours. 
































the scale height rather than the inclination. We therefore 
added a flare to the B model by lowering the scale height in 
the central rings. The resulting F model has a central scale 
height of 0.5 kpc out to a radius of 4.2 kpc. It then linearly 
increases to a maximum value of 1.6 kpc at a radius of 32.7 


kpc and stays constant after that. As can be seen in Fig. 27 


the addition of a flare improves the agreement with the data 
in the inner regions of the disk. 


On the other hand the F model still does not reproduce 
the forked shape of the observed emission in the outer parts 
of the disk. This forked shape is present on both the ap¬ 
proaching and the receding half, which seem to behave very 
similarly at first sight. In a first attempt to model this be¬ 
haviour we included a strong and abrupt line of sight warp 
in the F model. With an inclination of 79° in the outer part 
of the warp, the WFl model broadly matches the forked 
shape of the observed emission on the approaching half 
Indeed, although the model is slightly offset from the data 
in the very central channels and it does not include the faint 
tail that is observed above or on the back side of the disk, it 
is able to reproduce the general trend very well. On the re¬ 
ceding half this is not the case. If the position of the model 
matches that of the observed emission in the outer channels 
and the inclination is taken to reproduce the opening of the 
upper arm, the lower arm of the model opens up far too 
strongly. This is clearly visible in Fig. e.g. at 775 km 
s“^. A smaller line of sight warp could resolve this prob¬ 
lem, but also requires a change of the position angles in the 
warp to avoid now underestimating the opening of upper 
arm. Although this causes a clear offset between the model 
and the data at terminal velocities (at 898 km s“^ in Fig. 
[Z7] l we investigated this option with the WF2 model. The 
agreement of this model with the data in the lower arm is 
still very poor. Indeed, it now underestimates the opening 
of the data in the outer channels while still opening up a bit 
too strongly in the central channels. In fact, the only way to 
match the observed position of the lower arm in all reced¬ 
ing channels is with an almost edge-on orientation in the 
warp. This is demonstrated in Fig. |^by the WF3 model, 
with inclinations between 87° and 89° throughout most of 
the warp and dropping to 81° only in the outer tip of the 
disk (as shown in Fig.[29|. 

The two arms of the forked emission on the receding 
side thus seem to behave systematically different and can¬ 
not be modelled simultaneously with a simple, axisymmet- 
ric model. The lower arm behaves like an edge-on disk 
while the upper arm opens up and seems to suggest a sig¬ 
nificantly lower inclination. Unfortunately it is not possible 
to give a clear interpretation of this feature based on our 
data. Maybe only the back side of the receding half has a 
line of sight warp and the front side does not. On the other 
hand, the upper arm could also be caused by an extraplanar, 
arc-like feature related to the recent or ongoing interaction. 
Given the close companion and the stellar tidal streams that 
have been observed in the optical, it is indeed not unlikely 
that the gas disk of NGC 5907 is somewhat disturbed. As 
further modelling based on our data would degrade into 
guessing and speculating, we take the WF3 model as our 
(preliminary) final model and refer the reader to the future 
work of Yim et al. (in prep.). Their much deeper WSRT 
data will reveal the full extent of the atomic gas disk of 
NGC 5907 and provide additional constraints on its struc¬ 
ture. 


WF3: Warp app. + Edge-on rec. + Flare 



Fig. 28: Observed major axis position-velocity diagram of 
NGC 5907 (black and gray contours) overlaid with the fi¬ 
nal {WF3) model (green contours). Contour levels are -0.5, 
0.5 (1.5cr), 4.8, 9.6 and 14.4 mJy beam“^. The greyscale 
corresponds to the observations. 
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Fig. 30: Total Hi maps of the basic (B) and flaring (F) mod¬ 
els of NGC 973 as compared to the observed total Hi map. 
Contour levels are the same as in Fig. 


The variation of the radially dependent parameters of 
the model is shown in Fig. In contrast to the warped 
outer regions, the inner disk of NGC 5907 appears to be 
remarkably symmetric. The best fit values of the remaining 
parameters are listed in Table The inclination, position 
angle and scale height have errors of 1.0°, 1.5° and 0.3 kpc 
respectively. A full comparison of the final model to the 
observed data cube is given in Fig. A.6 of the appendix and 
the major axis position-velocity diagrams are compared in 


Fig. 28 
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Fig. 29: Variation of the radially dependent parameters of the final models. From top to bottom; the Hi surface density, 
the rotation velocity, the inclination, the position angle and the vertical scale height. Blue indicates the approaching side, 
red the receding side. For NGC4217, the parameters of the B+R2 model for the outer ring are also shown. The dashed 
vertical line shows (approximately) where the main disk meets the outer ring. For NGC 973, the flare from the F model 
(Section]^ is shown as the dashed red line. The uncertainties are represented by the error bars. 
























































































































Table 3: Radially constant parameters of the final models. 


Galaxy 

Center a 
(J2000) 

Center 5 
(J2000) 

^sys 

(kms"^) 

Velocity dispersion 
(kms"’^) 

Radial motions 
(km s-’^) 

NGC 973 

02''34"'20.F ± 0.4’’ 

-^32'’30'19.9" ±4.0" 

4770 ± 8 

8.2 ±2.0 

- 

UGC4277 

08‘'13"'57.0’’ ± 0.7” 

±52‘=38'53.9" ± 5.0" 

5370 ± 6 

10.4 ± 3.0 

- 

IC2531 

09''59"'55.5” ± 0.5” 

-29‘’37'02.4" ± 4.0" 

2455 ± 6 

11.2 ± 1.5 

- 

NGC 4217 

12''15"'51.1” ± 0.6” 

±49‘’05'29.4" ± 4.0" 

1022 ± 6 

19.1 ±3.0 

- 

NGC 5529 

14‘'15"'34.0” ± 1.0” 

±36°13'39.0" ± 8.0" 

2830 ±6 

12.3 ± 3.0 

-15 ±5 

NGC 5907 

15''15"'53.2” ± 0.9” 

±56°19'42" ± 15" 

664 ±7 

13.3 ±3.0 

- 


5. Constant scale height or flare 

As was discussed in Section our final models generally 
have a constant vertical scale height throughout the entire 
gas disk. A flaring Hi layer is only required for NGC 5907. 
Theoretically, however, it is expected that all spiral galax¬ 
ies have a flared atomic gas disk. Indeed, in the inner re¬ 
gions, the gravitational potential is dominated by the stellar 
disk. If the vertical velocity dispersion of the gas (and hence 
the internal pressure in the gas disk) is independent of ra¬ 


scale height of this model was only 0.1 kpc, which is far 
beyond our error bars (0.5 ± 0.2 kpc). A flare within these 
uncertainties is not in disagreement with the data, but does 
not significantly improve the model either (and was there¬ 
fore not discussed in Section 4.4.2| l. Hence, the arguments 
from the previous paragraph are also valid for IC2531 and 
NGC 4217. Within our uncertainties on the scale height, a 
flare can be added to the models without affecting the qual¬ 
ity of the fit, but its parameters are ill constrained. 


dius, the decreasing surface density of the stars with radius 

will result in an increasing scale height of the gas (e.g. |van| 6. Impact Of the Interactions 
1981|l. Beyond the stellar disk, the gravitational 


der Kruit 


potential is dominated by the dark matter and the flaring 
of the gas disk is directly related to the shape of the DM 
halo (e.g. |011ing||1995l|1996[|Becquaert & Combes||1997[ 


[O’Brien et al.|20i0d|b a|cHPeters et al.|2013|). 

We stress that these theoretical expectations are not in¬ 
consistent with our models. Indeed, as can be seen in Fig. 
29 the uncertainties on the scale height are generally sig¬ 
nificant and, in fact, leave room for a (moderate) flare. We 
illustrate this point for NGC 973, by replacing the constant 
scale height of 0.2 kpc in the final (B) model with a flare. 
In the resulting F model, the scale height linearly increases 
from 0.2 kpc in the center to 1 kpc at a radius of 45 kpc 
(dashed red line in Fig. 1^. From Figs. and it is 


clear that the B and F models are very similar. A flaring 
model is thus still fully consistent with the data and would, 
in fact, be more realistic from a theoretical point of view. 
The problem is that the details of the flare are generally ill 
constrained. Flares with different slopes or shapes would fit 
the data equally well and there is no way to tell which is 
best. For this reason we stick to a constant scale height in 
our final models, but note the reader that a flare is proba¬ 
bly present within the uncertainties on the (constant) scale 
height. The range of scale heights allowed by our uncertain¬ 
ties is indeed generally of the same order as the amplitude 


of the flares found in other galaxies, like NGC 4244 (Oiling 
[T9961 ) or the Milky Way ( jKalberla & Dedes|2008| ). 

Finally we remark that from Sections 4.3.2| and 4.4.2 
it might seem like we explicitely investigated the possibil¬ 
ity of a flare for IC2531 and NGC 4217 and discarded it. 
In fact, this is not entirely true. For both galaxies, only a 
certain type of flare was tested to model a specific feature 
in the data cube. In the case of IC2531, this flare did not 
significantly improve the model, but was not excluded by 
the data either. It was simply discarded to reduce the num¬ 
ber of free parameters. For NGC 4217, a flare was tested 
to reproduce the thinning of the emission in certain chan¬ 
nel maps. When this was achieved, the resulting model now 
appeared slightly too thin in other channels and was there¬ 
fore discarded. It should be noted, however, that the central 


Figure [^compares the radial variation of the main param¬ 
eters of the final model for each galaxy. The final param¬ 
eters that do not vary with radius are listed in Table and 
an overview of the main features of each model is given in 
Table Although there are, of course, differences between 
the individual galaxies, they all show normal behaviour and 
there are no clear outliers in the sample. 

As was discussed in Section the data cubes of 
NGC 4217, NGC 5529 and NGC 5907 contain a number of 
morphologic and kinematic indications of recent or ongo¬ 
ing interactions with satellite galaxies. From Fig. 29 it ap¬ 


pears that, apart from the radial inflow in the (outer) disk 
of NGC 5529, these interactions have not strongly affected 
the main Hi disks of the galaxies in question and their in¬ 
fluence is mostly limited to the position angles and inclina¬ 
tions in the outer regions. Indeed, the radial column density 
profiles and rotation curves of NGC 4217, NGC 5529 and 
NGC 5907 do not behave systematically differently than 
those of the rest of the sample and there are no strong asym¬ 
metries between the approaching and the receding halves. 
On the other hand, NGC 5529 and NGC 5907 do have a 
significantly stronger edge-on warp than the other galaxies. 
For NGC 4217 the edge-on warp has an unusual profile on 
the receding side, but its amplitude is similar to the warps 
of UGC4277 and IC2531. NGC 5907 also has a stronger 
line of sight warp than the other galaxies, but, on the other 
hand, NGC 4217 shows no line of sight warp in its Hi disk 
(only in the outer ring, as represented by the B+R2 model). 

In addition to these arguments, the galaxies in our sam¬ 
ple are all in agreement with the Hi mass-size relation that 
was found by Broeils & Rhee ( 1997| ). This can be seen in 
Fig. 32 and further suggests that the atomic gas disks of the 
‘interacting galaxies’ have not been strongly disrupted by 
the interactions. Following [Broeils & Rhee \\991) , we de¬ 
termined the Hi diameter Dhi from our final models as 
the diameter where the Hi column density falls below 1 
Mq pc“^. Since our models cannot adequately reproduce 
the outer ring of NGC 4217, we excluded this structure and 
only used the main disk of NGC 4217 to determine the data 
point in Fig. 32 We note that with the inclusion of the outer 














































Data B: Basic F: Flare 
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Fig. 31: Representative channel maps from the observed data cube of NGC973 (black and grey contours), the basic (B) 
model (green contours) and the flaring (F) model (red contours). Contour levels are -2.3, 2.3 (l.Str) and 6.2 mJy beam“^. 
The systemic velocity is 4770 ± 8 km s“^. 


ring, based on the B+R2 model, NGC 4217 is still in agree¬ 
ment with the relation from |Broeils & Rh'^ ( | 1997) 1. 

7. Conclusions 

In the frame of the HEROES project we present an anal¬ 
ysis of the atomic gas content of 6 edge-on spiral galax¬ 
ies. This work fits in a larger effort to derive the ra¬ 
dial and vertical distribution and the properties of the 
gas, dust, stars and dark matter in a sample of 7 massive 
edge-on spiral galaxies. New GMRT Hi 21-cm observa¬ 
tions of NGC 973 and UGC4277 are combined with re¬ 
reduced archival data of IC2531, NGC 4217, NGC 5529 


and NGC 5907. The atomic gas content of the seventh 
HEROES galaxy, NGC 4013, was recently studied accord¬ 
ing to the same strategy by |Zschaechner & Rand| ( |2015| l. 

The data cube of each galaxy is analysed through a set 
of detailed tilted-ring models. The primary goal of our mod¬ 
elling is to determine the main properties of the atomic gas 
disks such as the face-on column density distribution, the 
kinematics and the 3D geometry in an accurate and homo¬ 
geneous way. Eurthermore we identify and analyse peculiar 
features such as (potential) spiral arms and signs of recent 
interaction and separate them from the main gas disk. 

































Table 4: Main features of the final models 


Galaxy 

Edge-on warp 

LOS warp 

Flare 

Signs of spiral arms 

Signs of interaction 

NGC 973 

- 

- 

- 

- 

- 

UGC4277 

/ 

- 

- 

- 

- 

IC2531 

/ 

/ 

- 

Bright ridge in 
XV-diagram 

- 

NGC 4217 

/ 

- 

- 

- 

Coplanar ring outside 
main disk 

NGC 5529 

/ 

/ 


Bright ridge in 
XV-diagram 

Hi bridges to companions, 
extraplanar gas, 
radial motions 

NGC 5907 

/ 

/ 

/ 


Extended reservoir of 
low-density gas, 
asymmetric outer disk 



Log(D^,/ kpc) 


Fig. 32: The correlation between the total Hi mass and the 
diameter of the Hi disk for our sample (coloured symbols) 
and the sample of Broeils & Rhee ( 1997| l (grey dots, based 
on their Table 1). The black line shows the best fit linear 
relation from Broeils & Rhee ( 1997[l. 


Edge-on warps of varying strength are found in the 
outer regions of all galaxies except for NGC 973. In ad¬ 
dition we also identify moderate line of sight warps in 
IC 2531, NGC 5529 and NGC 5907. Evidence for a flar¬ 
ing gas disk is only found in NGC 5907, although, for 
the other galaxies, the uncertainties on the scale height 
leave room for a flare as well, in agreement with theo¬ 
retical expectations. 

The major axis position-velocity diagrams of IC2531 
and NGC 5529 show a bright Hi ridge. Our models fail 
to account for these ridges, although they adequately re¬ 
produce the other features in the XV-diagrams of these 
two galaxies. This suggests that the ridges are not part 
of an axisymmetric disk and supports the interpretation 
by [Kregel & van der Kruit ( 2004|l as prominent spiral 


arms. 


• Three galaxies in our sample show evidence of recent 
or ongoing interactions. We detect a (clumpy) coplanar 


ring of Hi gas just outside the main disk of NGC 4217. 
This structure has not been reported before and seems 
to be both kinematically and spatially offset from the 
main disk. We suggest this could be the relic of a re¬ 


cent minor merger event. Eor NGC 5907 Shang et al. 


(19981 already reported the detection of a compan¬ 
ion galaxy close to and a giant stellar tidal stream 
around the galaxy. Our modelling further reveals that 
the outer Hi disk is disrupted, with a strong (edge-on 
and line of sight) warp and large asymmetries between 
the approaching and the receding side. Einally, as al¬ 
ready noted by [Kregel et al. ( 2004| l, Hi bridges con¬ 
nect NGC 5529 to two nearby companions. Apart from 
these Hi bridges, we And that the gas disk is further 
surrounded (in projection) by a significant reservoir of 
extraplanar gas, mainly on the approaching side. A ra¬ 
dial inflow of 15 ± 5 km s“^ in the gas disk is possibly 
also related to the ongoing gas accretion. 

• The effect of these interactions on the gas disks of the 
galaxies in question is mostly limited to the orientation 
(inclination and position angle) of the outer parts. The 
column density profiles and rotation curves behave nor¬ 
mally and show no strong asymmetries between the ap¬ 
proaching and the receding halves. 

• Einally, a small companion galaxy is detected in the 
data cube of NGC 973, located at a projected distance 
of 36 kpc from the center of the main galaxy. We And 
no optical counterpart of this galaxy. 

In addition to the atomic gas, we are currently collecting 
CO line data to constrain the molecular gas content as well. 
An analysis of the dust and stellar content of the HEROES 
galaxies by fitting radiative transfer models to optical/NIR 
images is also underway. The next step in the HEROES 
project will be to At dynamical models to the rotation curves 
derived in this work to determine the distribution of the dark 
matter. Combining our analyses of all the different compo¬ 
nents, we will then ultimately investigate the correlations 
between them. The peculiar features identified in this pa¬ 
per will be an important factor in the interpretation of these 
results. 
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Appendix A: Comparison of the finai modeis 
to the observed data cubes 
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Right Ascension (J2000) 


Fig.A.l: Channel maps from the observed data cube of NGC973 (black and gray contours) and our final (B) model 
(green). Contour levels are -2.3, 2.3 (1.5 ct) and 6.2 mJy beam“^. The black cross indicates the center of the galaxy. The 
companion galaxy is indicated by arrows at 4531 and 4565 km s“^. 
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Right Ascension (J2000) 


Fig. A.2: Channel maps from the observed data cube of UGC 4277 (black and gray contours) and our final (B) model 
(green). Contour levels are -1.5, 1.5 (1.5cr) and 4.0 mJy beam“^. The black cross indicates the center of the galaxy. 
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Fig. A.3; Channel maps from the observed data cube of IC 2531 (black and gray contours) and our final ( W) model (green). 
Contour levels are -1.6, 1.6 (1.5cr) and 6.5 mJy beam“^. The black cross indicates the center of the galaxy. 
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Fig. A.4: Channel maps from the observed data cube of NGC4217 (black and gray contours) and our final (B) model of 
the main disk (green). Contour levels are -1.5, 1.5 (1.5cr) and 8.0 mJy beam“^. The black cross indicates the center of the 
galaxy. 









Declination (J2000) 



Right Ascension (J2000) 


Fig. A.5: Channel maps from the observed data cube of NGC 5529 (black and gray contours) and our final {W2R) model 
(green). Contour levels are -0.8, 0.8 (1.5cr) and 5.6 mJy beam“^. The black cross indicates the center of the galaxy. 
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Fig. A.6; Channel maps from the observed data cube of NGC 5907 (black and gray contours) and our final {WF3) model 
(green). Contour levels are -0.6, 0.6 (for better readability we take the lowest contour at 2a instead of 1.5 ct) and 6.4 mJy 
beam”^. The black cross indicates the center of the galaxy. 
































